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Determination of the characteristic of a 
non-linear resistor by harmonic excitation 


by prof. dr. C. J. Bouwkamp, Philips Research Laboratories, 


N.V. Philips’ Gloeilampentfabrieken, Eindhoven 


Consider a non-linear resistor with centre-symmetric current- 
voltage characteristic í = f(v) = —f(—v), and let this resistor 
be excited by a harmonic signal r = x cos(mwt), where the 
amplitude x covers a certain range (0 S x S a). Assume further 
that the amplitude g(x) of the first harmonie g(x}cos (wt) in 
the output current f(x cos (wt)) is measured in the same range. 
Then the question is: 
— Can the function f(x) be determined uniquely from the 
function g(x), where 0 S x S a; and if so, how can it be done? 
The answer to this question is of importance in relation to certain 
measuring techniques in experimental physics [I] and is re- 
ported in this note. 

If the resistor is excited by rv = x cos p. where p = wt, the 
output is f(x cos @), and the amount of the first harmonic in 
this output clearly is 
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In virtue of the symmetry of the cosine function and f(v) being 
an odd function of v, we can replace the upper limit by n/2 
and multiply the expression by 4. Therefore, the two functions 
fand g are interconnected by the following integral equation: 


nii 

[ f(x cos p) cos rPdp= jl e(x), Osx<a 

Li 4 

This integral equation is closely related to Abels integral 
equation, as can be seen by introducing the new variable of 
integration py = X COS @. 
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On substitution we get 


LOON Wzn 


0 ve k 


The solution of this latter equation is known. Under certain 
mild conditions as to smoothness of f and g. not to be specified 
here, the integral equation has the unique solution 


F 
d ( g(x) x’ dx 


Ì 
Le == EE de 
y tú vr-x 


and this answers our question. 


(0O< y <a) 


This solution in closed form though is not very attractive to 
the practical engineer if g is not known analytically. An alter- 
native solution by means of power series or the like is more 
practical, especially so if the resistor is slightly non-linear. 
Therefore, let us assume that 


glo)=se, gent 


where the right-hand expression may be an infinite convergent 
series, or a least-square polynomial approximation to measured 
values of g, or a Chebyshev approximation, etc. If this is 
substituted in the closed-form solution, the result is an ex- 
pression of the same type, namely 


LD) - y gen nl (n +1)! é p= +1 
(2n + 1)! 


where the summation is extended over the same set of values n 
as in the expression assumed for g. In the case of infinite series. 
the two power series have the same radius of convergence. 

More explicitly, if the measured values of g can be represented 
by a polynomial approxtmation, say 


90) Zeg Herk Hex +eax! 
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then the required characteristic is of the form 


M= a Ke nd Ban Ot nl 


of course approximately, and over the same range. 


Finally, it should be noted that in practice we need not invoke 
the theory of Abels integral equation at all, apart perhaps from 
the question of uniqueness. In fact, if we take f (rv) = v” (vr = 0). 
where m is not necessarily an integer, then the corresponding 
glx), x 2 0, is apparently given by 


if m 
4 nij ple) 
Xx g(x) =— | gt P - en ON nn 
rrd 
À À ) 
and therefore. if 
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then 


The author is indebted to Professor Butterweck, Technological 
University Eindhoven, for drawing his attention to this problem. 
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strength characteristic of n-type gallium arsenide using various high- 
power microwave techniques. Trans. LE.E.E., ED-14, 505-511, 1967. 


On the stability of a circuit with negative resistance with 


non-linear characteristic 


by prof. ir. J. W. Alexander, peift University of Technology, Dept. of Electronical Engineering 


Summary: For an amplifier with a negative resistance the operating point on the curved characteristic of the element 
has to be stable. A criterion for the stable points on the curved characteristic is given, without the addition of an extra 
capacitance or inductance. This criterion runs as follows: dE/d/ must be positive, where E and / are the adjustment 
voltage and ditto current respectively. When the circuit contains also other elements than only resistances the amplifier 
circuit is stable when in addition to the adjustment stabilization the poles must have a negative real part. 


1. Introduction 


Á negative resistance always has a curved characteristic with a 
negative slope on a part of the curve. For operation as an 
amplifier we have to choose our operating at any point on the 
negative slope of the characteristic. For stable operation this 
operating point has to be a stable one. In the discussion on this 
problem one usually takes into account the existence of an 
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extra element, either an inductance or a capacitance. In the 
following an other method is used, without the addition of an 
extra element. 


2. Criterion of stability 


We will start from the fact that a circuit consisting of a d.c. 
source with voltage E and a positive resistance R is always 
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Fig. 1. Circuit with a positive and a linear negative resistance. 


ke 


Fig. 2. Characteristic of a non-linear negative resistance with different 
laad lines. 


stable. That means that a small change dE in the voltage E 
gives a small change d/ in the current, in such a way that dE/d/ = 
+ R, a positive value. 

For a linear negative resistance — R we should have E = — IR, 
sodE/df < 0. This is a negation of the stable situation dE/d/ > 0, 
consequently an unstable situation, in agreement with the im- 
possible situation E = — IR. 

When the circuit consists of a positive and a negative resistance 
in series, according to Fig. 1, then dE/d/ = Rr. 

When R > r, dE/dl is positive and the circuit is stable, when 
R < r however, it is unstable. 


3. Circuit with a non-linear negative resistance 


Á negative resistance always has a curved characteristic with a 
negative slope on a part of the curve. Let Fig. 2 represent the 
U vs. F-diagram of the non-linear element. 

In this figure A is the operating point chosen in the inflexion 
point of the characteristic, where 


ER, lo = Uo ( 


where E, is the adjustment voltage and tan(«,) = R, 
In the figure another situation is also shown, where 


E‚ Rl, = U) with tan («,) = R, 

With 

U= U, +uand/= 1, + î (2) 
we may describe the curved characteristic by the equation 
u=-ri+ bi (3) 


from A till beyond the extreme values. Here r and h are positive 
and tan(#) = r. We see in the figure that R‚ > rand R‚ <r. 

The points of intersection of the characteristic and the load 
line follow in the first situation from 


u=-—ri + bi together with u = — Ri. 
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This results in í = O and —R‚ = —r + bi,orbi =r-R,. 
Now R‚ > r, so we have only one intersection point A. 
For the second case we find in the same way 
i= Oand bi =r—R, 
or the three intersection points A, B and C. 


For the discussion of the stability we proceed as follows: 
In the first situation Fig. 1 shows that 


EE, =RI+ U (4) 
Substituting (1), (2) and (3) in (4) we obtain 
E — IoR, — Up =(R‚ — mi + bi. 


For the adjustment stability we get 
dE/dI = dE,‚/di = R‚ —r + 3 hi 


for F, and U, are constant, 

The only intersection point is i = 0. As R‚ > r the point A 
is stable. 

But for the second case with R‚ < r we get 


dE/dI = dE‚/di = R‚ —r + 3 hi. 

For A: i = Oand R‚ < r, consequently A is unstable. 

For Band C: bi? = r — R,, so 

dE/dI = R‚ —r + 3(r — R‚) = 2(r — Ro). 

As R,‚ < r, B and C are stable points. 

Consequently: For operating in the inflection point A the 
positive resistance R must be greater than the negative resist- 
ance r. 


4. Situation with a negative conductance 


Till now we have studied the element with the negative resist- 
ance, according to Fig. 2, or the element with the current- 
controlled characteristic. 

The other element with the voltage-controlled characteristic, 
Le. the same form as in Fig. 2, but with U and / interchanged, 
has the criterion of adjustment stability: dA/dU” > 0, where A 
is the adjustment current. This results in the similar condition 
that for-operating in the inflection point A the positive con- 
ductance G must be greater than the negative conductance g 
of the element. 


5. Amplification with a negative resistance and other elements 


In this case the requirement for the adjustment stability is not 
the only one which has to be satisfied. The poles of the charac- 
teristic equation of the circuit must have a negative real part, 
or the Hurwitz criteria are to be satisfied also. 

To demonstrate this we examine the circuit, shown in Fig. 3, 


Fig. 3. Circuit for amplification with a current-controlled negative 
resistance. 
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for amplification with a current-controlled negative resistance. 
For the adjustment stability we get from / = El(r, + rn, —F) 
that 


ptn >r (5) 


is required. 
For the time-depending e we get from 


u=er (l —rpO/ir, + r, — r — (r, + 7) rpC} 
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the characteristic equation: 

he Pen Elft APC =D. 

For a stable situation it is well-known that the pole must be 
negative, or 

A (6) 


Now (5) and (6) cannot be realized at the same time, conse- 
quently: the circuit according to Fig. 3 with a current-controlled 
negative resistance is not possible as an amplifier. 


Niet-lineariteit van de lawine-looptijd-oscillator 


door dr. ir. K. Mouthaan, Natuurkundig Laboratorium, 
N.V. Philips’ Gloeilampentabrieken, Eindhoven 


Summary: Non-linearity of the Avalanche Transit-Time Oscillator. 


High-power solid-state generation of microwaves has become possible with the Avalanche Transit-Time Oscillator. 
The oscillator is based on the dynamic negative-resistance properties of semiconductor diodes that are reverse- 
biased into avalanche breakdown. Optimum design conditions are sought by analysis of the non-linear operating 
characteristics. Simultaneous operation at two harmonically related frequencies is shown to permit well-defined 
design of the oscillator, and to yield high output power and low frequency-modulation noise. Theoretical tuning 
conditions and curves for output power as a function of the bias current are presented. 


De lawine-looptijd-oscillator is een recent tot ontwikkeling ge- 
komen vaste-stof microgolfbron. Frequenties van 1 GHz tot 
meer dan 300 GHz en nuttige vermogens tot enige watts in 
continu bedrijf zijn experimenteel gerealiseerd. 

Voorgesteld door Read [1], berust de oscillator op looptijd- 
negatieve weerstand van een in lawinedoorslag bedreven half- 
geleiderdiode. Een succesvolle diodestructuur is die met asym- 
metrische abrupte p-n-overgang. Het profiel van het elektrische 
veld binnen een dergelijke in sperrichting aangesloten structuur 
is schematisch voorgesteld in fig. 1. Het verontreinigingsprofiel 
is zo dat bij voldoend hoge sperspanning enerzijds lawine- 
doorslag optreedt op de p *-n-overgang ten gevolge van stoot- 


le Ld 


Fig. 1. Schematische voorstelling van het profiel van het elektrische 
veld in een lawine-looptijd-diode met asymmetrische abrupte p-n- 
overgang. In het gearceerd aangegeven doorslaggebied ter lengte /, 
is de veldsterkte van de orde van 200... 500 kV/cm (silicium). Tussen 
opeenvolgende botsingen met het kristalrooster kunnen vrije lading- 
dragers in dit gebied voldoende energie uit het veld opnemen om 
stootionisatie te bewerkstelligen. In het driftgebied ter lengte /, blijft 
de veldsterkte groter dan 10... 20 kV/cem. 
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ionisatie; anderzijds strekt de depletielaag zich dan uit tot aan 
de n"-laag. In de gehele depletielaag is daarbij de veldsterkte 
voldoende hoog om de in de doorslag gegenereerde elektronen 
de n-laag met veldverzadigde driftsnelheid te doen doorlopen. 
De tevens gegenereerde gaten verdwijnen direct in de p*-laag. 

Een boven de gelijkspanning aanwezige wisselspanning modu- 
leert de in de doorslag gegenereerde hoeveelheid elektronen. 
In het als driftgebied fungerende gedeelte van de n-laag onder- 
gaat ten gevolge van looptijd de gemoduleerde elektronen- 
stroom een continue fasedraaiing ten opzichte van het wissel- 
veld. Voor voldoend grote looptijd resulteert daaruit een nega- 
tieve weerstand. Heeft de modulatie van de elektronenstroom 
een tijdvertraging gelijk aan een kwart periode ten opzichte 
van het wisselveld, dan treedt negatieve weerstand op voor elke 
waarde van de looptijd. Dit geval van inductieve stroomgene- 
ratie is karakteristiek voor de lawine-looptijd-diode. 

Is V‚(#) de boven de gelijkspanning aanwezige wisselspanning 
over het doorslaggebied, en is daar «(E) de van het elektrische 
veld afhankelijke ionisatiesnelheid, d.i. het aantal gat-elektron- 
paren dat één ladingdrager per eenheid van lengte doet ont- 
staan, dan geldt voor de in het driftgebied geïnjecteerde con- 
vectiestroom J‚(f) bij benadering: 


br | 
J.W) =J, exp en | V_ (@) el (I) 
T 


a 
C 


Hier is «/ (>0) de bij het aangelegde de-veld geëvalueerde 
afgeleide van « naar het elektrische veld en z, is de tijd benodigd 
voor het doorlopen van het doorslaggebied. De de-instel- 
stroom Jo = <SdJ,lt)> wijkt voor willekeurig signaalniveau 
V‚(t) af van de stroom J,; wordt de diode gevoed uit een 
constante-stroombron, dan is J, een onafhankelijk variabele. 
Uit vel. (1) blijkt het inductieve karakter van de stroomgene- 
ratie; voor kleine waarden van V(t) kan dit karakter worden 
voorgesteld door een zelfinductie L, = t‚/(2x,J). 
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De niet-lineariteit van de stroomgeneratie volgens vgl. (1) 
compliceert het gedrag van de lawine-looptijd-oscillator in 
belangrijke mate. Ten behoeve van het inzicht in de werking 
en het verkrijgen van ontwerpcriteria wordt in het volgende een 
analyse beschouwd van het niet-lineaire gedrag van de diode 
plus aangekoppeld mierogolfcircuit. De vanwege het aange- 
legde de-veld in het driftgebied geïnjecteerde convectiestroom 
Jt) gedraagt zich daar als een niet-dispersieve golf: 


Jalx, ĳ) = Jal0. t—x[va) 5 Talt xv) 


Hier is J‚,(x,t) de convectiestroom in het driftgebied en v, de 
veldverzadigde driftsnelheid van de ladingdragers. De wissel- 
spanning over het driftgebied zij V‚(4); de totale wisselspanning 
over de diode is dan V‚(t) + V‚lt). Behalve de convectiestroom 
is ook de verschuivingsstroom continu op de overgang van 
doorslaggebied en driftgebied; in het eendimensionaal onder- 
stelde model is de som van convectiestroom en verschuivings- 
stroom in het driftgebied divergentievrij. 

Voor cyclische werking van de oscillator kunnen alle tijd- 
afhankelijke grootheden worden voorgesteld door een fourier- 
reeks. Zo schrijven we voor J‚‚(4): 


fbd | +5 Bret" 4 cl @) 
en voor V‚(t): 
PbN Wie" ec) 3) 


Hier is i de imaginaire eenheid en de notatie c.c. duidt op de 
complex geconjugeerde term. De stroomamplitudes f,, worden 
bepaald door substitutie van vel. (3) in vgl. (1); ieder van de 
stroomamplitudes f,, hangt dan (niet-lineair) af van alle span- 
ningen V,…. Van belang is dat de stroom J‚(f) positief definiet 
js; daaruit volgt dat |#‚| <1. Is nu C, de capaciteit van het 
doorslaggebied, dan wordt de mde harmonische component 
van de totale stroom op de overgang van doorslaggebied en 
driftgebied gegeven door 

Form = Êm Teo + iMwC, Vi (4) 
De impedantie van de diode voor de frequentie mo, Zaioae, m — 
(Win + Var) me kan dan op grond van het aangenomen 
model worden berekend als - 


i | 
nw Ì + B den immwC, Kas 


Laiode, Er _ nr 


| 5) 
l(L— et"°)/ immo 


Ee ‘ = en 
imwC, | Li imaC, Vn |B Jo 


Hier is C, de capaciteit van het driftgebied, C is de totale capa- 
citeit van de diode volgens C = C,C‚/(C, + C), en 0 is de 
loophoek voor de frequentie w volgens 8 = wl;/v,. Het looptijd- 
effect wordt uitgedrukt door de tweede term in vgl. (5). Behalve 
door de geometrische parameters C, C‚ en @, wordt het im- 
pedantiegedrag van de diode geheel bepaald door de verhou- 
dingen imoC,V/P Je, dus door de verhouding van ver- 
schuivingsstroom en convectiestroom in het generatiegebied. 
Is de generatie op de frequentie » inductief, dan geldt f, = 
iV,|B‚/V‚,| en vertoont Zaioge,, een parallelresonantie voor 
de frequentie waar |wC,V‚‚l = |P‚Jool- De diode-impedantie 
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Ziioge, iS dan in principe beneden de resonantiefrequentie in- 
ductief en passief, boven de resonantiefrequentie capacitief en 
actief. 

De werking van de oscillator wordt bepaald door de voor- 
waarde Zisoae m + Zem = 0 voor alle m, waarin Z, de impe- 
dantie van het lineair en passief onderstelde externe circuit voor 
de frequentie mw voorstelt. De diode-impedanties Zasae, m 
hangen af van de instelstroom en het signaalniveau, terwijl de 
impedanties Z,„, niet afhangen van deze grootheden. Een zelf- 
startende oscillator voor bedrijf op een vooraf gespecificeerde 
frequentie en instelstroom vereist daarom in het algemeen een 
complex, zo niet onpraktisch, microgolfcircuit. Een bijzondere 
mogelijkheid bestaat nu indien voor waarden van J., boven de 
start-oscillatiestroom de spanningen V,, continu aanliggende 
waarden kunnen aannemen zodanig dat de verhoudingen 
moC Vn /B Jo constant blijven. In dat geval blijven de diode- 
impedanties Zia. „ en dus de oscillatiefrequentie constant bij 
variërende J,, en wordt een welgedefinieerd ontwerp van het 
aan de diode te koppelen microgolfcircuit mogelijk. Van bij- 
komend belang is dat dan tevens parasitaire frequentiemodu- 
latie ten gevolge van fluctuaties in J vermeden wordt. 

Om de genoemde mogelijkheid te onderzoeken, schrijven we 
de vergelijking Ziege m + Zem — 0 met voordeel als een ver- 
gelijking voor de verhouding imoC, Val Pd co: 


moC VenlBeso + Pm =O (6) 
Hier is: 

pn = Ben + Zand Lem + WimmC) 

en 


Zin = (fimoC,) {1 (1 —e*"”/im0) 


Zoals reeds opgemerkt is de stroomgeneratie voor voldoend 
lage waarden van V‚(f) inductief, dus B, = —iV,,lP/V.ol- 
Oscillatie is voor lage signaalniveaus derhalve slechts mogelijk 
indien p‚ reëel is, hetgeen de afstemvoorwaarde 


l 


wC 


C 0/2 
Im Zi= + (R,, +Re zjn u. Ed (7) 


C, sin(4/2) 


oplevert. Hierin is Im Z,, het reactieve deel van Zen A, + 
Re Z,, het resistieve deel van Z,, opgesplitst in een parasitair 
deel R,‚ en een nuttig deel Re Z,,. Voor reële waarden van p, 
geldt | 


p‚= 1 + Q, sin“(0/2)/(0/2), (8) 


waarin Q, = {WC,(R, + Re Z‚)} * een maat is voor de 
kwaliteit van het oscillatorcircuit op de frequentie w. Daar voor 
een passief extern circuit p,‚ > 1, zien we dat aan vgl. (6) slechts 
voldaan kan worden voor frequenties boven de resonantie- 
frequentie, terwijl voor gegeven Q, de instelstroom een vol- 
doend hoge waarde moet hebben om de oscillatie inderdaad 
te doen starten. 

Een oplossing als hierboven voorgesteld blijkt nu voor een 
willekeurig signaalniveau niet mogelijk te zijn indien de stroom- 
amplitudes f,, in essentie slechts bepaald worden door de 
grondgolf-spanning V,,. Naar orde van analytische en ook in 
de praktijk realiseerbare gecomp.iceerdheid, onderzoeken we 
derhalve de mogelijkheid waarbij de stroomamplitudes #,, in 
essentie bepaald worden door spanningen V,, op twee Îre- 
quenties, w en nw. De spanningen VW, voor k + ln zullen 
daartoe klein blijven ten opzichte van V,, en V,,; volgens 
vgl. (6) vereist dit dat de coëfficiënten p,‚ voor k # ln klein 
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gemaakt worden ten opzichte van p, en p „„ Met de spannings- 


normalisatie £, = Zia, V,/mot, worden dan de relevante 
oscillatorvergelijkingen: 
PADS Jep, valde 
BD, = (r?p‚/p) Jao, tart ol Teo) (9) 
waarin Jo, arro =S (CA/2a, vj) (W°/p,) die waarde van J,, voor- 


stelt waarvoor de klein-signaal-impedantie van het oscillator- 
circuit voor de frequentie w verdwijnt. De grootheid A is de 
oppervlakte van de dwarsdoorsnede van de diode. 

De stroomamplitudes f, worden dan voorts bepaald door de 
fourierontwikkeling van 


exp{ —(D,°°' + B*ele)} x exp{ — (Dele! + D*erirot), 


Daar slechts relatieve fasen van belang zijn stellen we 7, posi- 
tief reëel; de complexe spanning £, schrijven we als produkt 
van grootte en fasefactor: 7, — v‚e'**, De gevraagde fourier- 
ontwikkeling vinden we als een som van mengprodukten van 
de door de eerste exponentiële functie gegenereerde harmoni- 
schen van de frequentie w en de door de tweede exponentiële 
functie gegenereerde harmonischen van de frequentie nw. We 
vinden dan: 

"00 

ae} Í Ireen 1, 2v 5 Zen We (2) elP en 
B a ad eid (10) 

y (DPP TL (@o) 1, (2) e°v= 


De functies /, (x) zijn de gemodificeerde Besselfuncties van de 
eerste soort, de orde v en het argument x. Onmiddellijk volgt 
uit vgl. (10) dat f, voor willekeurig signaalniveau slechts een 
constante fase ten opzichte van v, kan behouden indien òf 
W, ne 0, òf mn 

In dat geval wordt tevens ff, reëel en vereist oplossing van 
vgl. (9) een reële waarde van p,. Deze eis levert voor de fre- 
quentie nw een afstemvoorwaarde-overeen kqmstig die voor de 
frequentie w, dus vgl. (7) voor de frequentie ne. De bijbehorende 
reële waarde van p, wordt gevonden als vgl. ($) voor de fre- 
quentie nw. 

Een belangrijke te bepalen grootheid is het door de oscillator 
op de frequentie w geleverde nuttige vermogen. Door met 
vgl. (6) de totale stroom volgens vgl. (4) uit te drukken in 
BJ cos Vinden we dat voor reële p‚: 


ES 
P, E e - Jo Q, oe BB Re Za (LI) 
2 012 
Slechts de stroomamplitude fl, is hierin de nog nader te bepalen 
grootheid, Daar |ff,| S 1, levert vgl. (11) direct het maximaal 
haalbare uitgangsvermogen. De optimale loophoek wordt ge- 
geven door } = 0,74. Opgemerkt zij dat het maximaal haal- 
bare uitgangsvermogen niet slechts gelimiteerd wordt door de 
eindige haalbare waarde van het produkt Q,* Re Z,,, maar 
tevens en vooral ook door de eindige toelaatbare waarde van 
de hoogfrequent spanningszwaai over de diode. In fig. 2 is 
uitgezet het door numerieke oplossing van veln. (9) en (10) 
verkenden relatieve uitgangsvermogen p‚… = (BJolJao san) 
= v‚° als functie van J olden san » met p/p, als parameter voor 
een bweelmqnenticentbinatie met n = 2. Als opvallende eigen- 
schap zien we in fig. 2 voor iedere waarde van p/p, het bestaan 
van twee modi van oscillatie: een modus met hoog vermogen 
en een modus met laag vermogen. Deze twee modi van oscillatie 
komen respectievelijk overeen met , = zen wy, = 0, hetgeen 
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Fig. 2. Relatieve waarden van het op de frequentie w geleverde nuttige 
vermogen als functie van de relatieve instelstroom bij gelijktijdig be- 
drijf op de frequenties w en 2w. De parameter P‚/p, is een maat voor 
de relatieve kwaliteit van het oscillatorcircuit op de frequenties 
en Zw. 


zich laat verklaren door het feit dat de oscillaties op de fre- 
quenties w en 2, hoewel gekoppeld, ten opzichte van elkaar 
nog twee verschillende fasen kunnen hebben. De start-oscillatie- 
instelstroom kan verschillen van J … „ en verschilt voor de 
twee modi van oscillatie, behalve indien p‚/p, = 4. Ingeval 
P/P, = 4 komen de twee vermogenscurven voor p‚… — 0 
samen bij JoJo aan = |; de modus met hoog vermogen kan 
dan slechts bereikt worden via een sprong in oscillatie-ampli- 
tude. De modus met hoog vermogen heeft de laagste waarde 
van start-oscillatie-instelstroom, hetgeen praktische realisatie 
van deze modus vereenvoudigt. Het hoogste vermogen wordt 
voor beide modi verkregen voor de grootste waarden van PalPy; 
daartoe dient de kwaliteitsfactor Q, van het oscillatorcircuit 
voor de frequentie 2w hoog gekozen te worden. Opgemerkt zij 
dat tweefrequentiebedrijf voor waarden van n= 2 minder effi- 
ciënt blijkt dan bedrijf met n = 2. Voorts vormt een uiteindelijk 
oscillator-ontwerp nog een compromis voor wat betreft toelaat- 
bare hoogfrequent spanningszwaai en voor wat betreft ther- 
misch toelaatbare waarde van de dc instelstroom. 

Voor een meer volledige beschrijving van de hier gegeven 
theorie en haar consequenties zij verwezen naar [2] en [3], waar 
tevens verdere verwijzingen naar de theorie van lawine-looptijd- 
oscillatoren gevonden kunnen worden. Opgemerkt zij nog dat 
met in het Natuurkundig Laboratorium in de groep van dr. 
D. de Nobel gemaakte dioden en constructie van de oscillator 
volgens de hier gegeven theorie, een nuttig vermogen van 1,75 W 
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in continu bedrijf op 5 GHz met een conversie-rendement van 
99 kon worden bereikt. 


a | 
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De toepassing van drie gekoppelde transmissielijnen 


door ir. R. Blommendaal, Christiaan Huyaenslaboratorium N.V, Noordwijk 
Dedicated to Protessor J. P, Schouten on the occasion ot his rebrement, 


Summary: The application of three coupled transmission lines. 


The theory of three coupled waveguides is outlined and applied to three waveguide rings bent in the E-plane. The 
input and output of the central ring are connected. This configuration can be used as a rotary joint if a gap is cut : 
at the centre line of the middle ring. Expressions are given for the transmission as a function of rotation and fre- Nes 


quency. 


|. Inleiding 


De theorie van de richtkoppeling tussen twee transmissie- 
lijnen is zeer goed bekend. In de literatuur is een overvloed aan 
stof te vinden over koppeling van golfgeleiders, coaxiale lijnen 
en ‘stripline’, waarbij zowel de continue als de discrete koppeling 
wordt behandeld. De golven in de secundaire lijn hebben steeds 
een faseverschil van 90° met die in de hoofdlijn en de ampli- 
tudes kunnen gevonden worden als functie van het koppel- 
mechanisme en de lengte van de gekoppelde lijnen. 

Het algemene geval van gelijktijdige en continue koppeling 
tussen een willekeurig aantal transmissielijnen is ook behandeld, 
maar slechts weinig bekend. Enerzijds komt dit omdat de op- 
lossingen snel ingewikkeld worden en niet kunnen worden af- 
geleid uit de theorie van twee transmissielijnen, anderzijds om- 
dat er weinig aantrekkelijke toepassingen bekend zijn die werk 
op dit gebied zouden kunnen stimuleren. 


De theorie van drie gekoppelde golfgeleiders zal beknopt 
worden beschreven en worden toegepast op een geval waarbij 
de golfgeleiders in het E-vlak tot ringen gebogen zijn. Ingang 
en uitgang van de middelste golfpijp worden daarbij met elkaar 
verbonden. Deze configuratie is met een kleine wijziging bruik- 
baar als draaikoppeling indien de middelste ring op de neutrale 
lijn gedeeld wordt. Uitdrukkingen zullen worden gegeven voor 
de transmissie als functie van de frequentie en de stand van 
in- en uitgang. 


E‚ (x) 
6 3 
| E Ge _ . Er ed X I Je TSE + Ì 
Ì 3 3 
È 5) ede Jed: +: p) iser Í 
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2. Drie gekoppelde lijnen 


In fig. 1 zijn drie gekoppelde lijnen geschetst. Gedacht wordt 
aan rechthoekige golfgeleiders met golfvoortplantingen in de 
TE, ‚-mode, die via de smalle zijde gekoppeld zijn. Voorts wordt 
aangenomen dat de koppeling continu verdeeld is, de golven 
zich alleen in de positieve x-richting voortplanten en dat de 
tijdafhankelijkheid exp(jwt) is. De volgende vergelijkingen kun- 
nen dan worden opgesteld: 


dE mas 

Ed iPE, Je (E, 0 E.) 

dx 

dE. PR | 
S= jE JE LE, + 2E, 4E) (1) 
dx î 

dE 


) 


SPE — ie lEr + B 
dx 

waarin: f — voortplantingsconstante 
c = koppelfactor 


Ì 


Oplossingen van de vorm E‚= exp(—jrx) kunnen worden ge- 
zocht en blijken te zijn: r, = f‚r, = f + cenr;z =P + 3c. 


Met behulp van de randvoorwaarden: E,(0) = 1, E‚(0) = 0, 
E,(0) = O en £,(0) = 0, E,(0) = len E,(0) = 0 kunnen de 
integratieconstanten worden bepaald en daarmede de coëffi- 
ciënten van de strooimatrix: 


me E‚ (0) 
= E‚ (0) (2) 
ACH LE A 
se te E, (0) 
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Fig. 1. Drie gekoppelde lijnen. 
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Fig. 2a, b. Amplitude en fase als functie van ex; 
E‚ (0) = 1, E, (0) = E, (0) =— 0. 


Omdat gerichte koppeling is verondersteld, is hier volstaan met 
een 3 x 3-matrix. Het verloop van amplitude en fase als functie 
van cx is gegeven in fig. 2a resp. 2b voor het geval E,(0) = 1 en 
E,(0) = E,(0) = 0. Voor geen enkele waarde van ex kan vol- 
ledige vermogensoverdracht naar de lijnen 2 en 3 worden be- 
reikt. 


3. De ringvormige configuratie 


De gevonden strooimatrix kan worden toegepast op de con- 
figuratie van fig. 3. Alle lijnen zijn gebogen met straal r‚ terwijl 
de middenring doorloopt. Het geheel kan worden opgevat als 
twee zespoorten ter lengte van rb, en rÖ,. Indien alleen poort | 
met amplitude = 1 wordt gevoed en de uittredende golf in 
poort n wordt aangeduid met b,, kunnen met behulp van (2) 
de volgende vergelijkingen worden opgesteld: 


Da X, B, Ni l 
bs =| ò, P, bu (3) 
be “B a bz 
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bo a, Ps, 72 B, 
b, e= B, Ö, P, bs (4) 
ba Ya P» 3 Ô 
met 
P‚+3p,+2 pp, tl, 
6 6 (n=lof 2) 
il Ip 41 
B R, TLG 
3 3 (5) 
p ederd, p eier û‚) 
es 2 
(0 <= Ö, ZS 2r) 
R, — edfr6: É R, =g Pr) 
Uit (3) en (4) volgt direct: 
B, Ba + rra Db + (Bô, + 1 B2) bs es | 
(Ö, Ba + B, 7) ba + (O0, +,B,-Ib; =-f, (6) 


B, Ba + a, va) ba + (B, 0, Haf) bs-b, =—y, | 


Hieruit laat zich b, oplossen als functie van B, 8, en c. In dit 
korte bestek zou het te ver voeren dit algemene geval uit te 
werken. Beperken wij ons eerst tot het geval dat 9, = 2r, dan 
geldt «, = ò, = len fl, = 7, = 0, zodat voor volledige ver- 
mogensoverdracht van | naar 6 moet gelden: 


= + ef (7) 
l_ó, 
Aan deze voorwaarde kan worden voldaan indien: 
WIR, + IpiR, pi + pk, — p-2 = (8) 
WIR, + R‚-3#0 | 


Mogelijke oplossingen: 


B So Med (9) 
en 


R‚=1 Bel 


Pi = =l 
R‚ = —j 
De laatste oplossing kan worden gerealiseerd met de kortste 


lijnlengte en komt bovendien overeen met de voorwaarde voor 
volledige overkoppeling bij twee lijnen. Beschouwen we dit 


Fig. 3. De ringvormige configuratie. 
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Fig. 4. De transmissie in amplitude en fase als functie van de stand 6. 


geval, waarbij cr = l/4en fr = Mm + 1/4 dan volgt uit (6) voor 
willekeurige 0, = Ö: 
be | ee eÎ Pi dam 


bm 


Dieos Aisne 4Ijsin 6-2 
4 d 2 


(10) 


waarin m geheel voor de hartfrequentie. 


4. De toepassing 


Het verloop van de transmissie van poort | naar poort 6 als 
functie van 0 is gegeven in fig. 4 (m = 5). De amplitude is 
constant binnen + 0,25 dB, terwijl de fase maximaal 17” van 
een lineaire toeneming afwijkt. Deze uitvoering is geschikt als 
draaikoppeling indien de middenring op de neutrale lijn gedeeld 
wordt zoals in fig. 5. Beide helften kunnen dan vrij gedraaid 
worden ten opzichte van elkaar. De variatie van de amplitude 
ten gevolge van het ronddraaien kan nog verminderd worden 
indien men de voortplantingsconstante van de middenring 
wijzigt. Dit kan gedaan worden door de breedte van de golfpijp 
te veranderen. wat door de aanwezigheid van deling eenvoudig 
is. De theorie van dit geval is ook te vinden op de bovenbe- 
schreven wijze, maar is zeer bewerkelijk. 

De ringvormige draaikoppeling heeft het voordeel dat zij op 
eenvoudige wijze uitgebreid kan worden tot meerdere kanalen, 
door identieke exemplaren te ‘stapelen’. Bij de gebruikelijke 
coaxiale draaikoppeling is dit niet mogelijk, al kan men een 
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Fig. 5. Toepassing als draaikoppeling in rechthoekige golfpijp. 


systeem met twee kanalen realiseren door de binnengeleider 
ook als coaxiale lijn uit te voeren. Het ruime gat in de ringvormige 
draaikoppeling maakt het bijv. mogelijk alle naar boven lopende 
lijnen door de gaten te voeren, en de naar beneden lopende aan 
de buitenzijde te houden. Naast het nadeel van verandering 
van de fase met de stand heeft dit soort draaikoppelingen een 
beperkte bandbreedte: 5 à 10%. Het gedrag als functie van de 
frequentie kan eveneens uit de boven gegeven uitdrukkingen 
worden berekend. Er treden in de 0’-stand stopbanden op nabij 
frequenties waar de weglengte door de middenring een oneven 
aantal halve golflengten is. Een praktische uitvoering voor 
frequenties van 8800 …. 9200 MHz heeft de bruikbaarheid 
van de rekenwijze aangetoond. 
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Modern and classical antenna design techniques 


by dr. ir. M. S. Afifi, Deltt University of Technology, now with Directorate for Research and 
Engineering of Page Communications Engineers, Inc, Washington D.C. 


Summary: A short review is given of the history of antenna design techniques. Some viewpoints are explained for 
modern techniques in comparison with the classical ones. More attention is given to frequency-independent an- 
tennas and the necessity of finding simple mathematical concepts in order to comply with the requirements of 
evolving physical pictures. 


1. Introduction are formulated in a compacted form by Maxwells equations. 


The antenna is a device for generating these waves with given 
Electromagnetic waves are controlled by physical rules which specifications which are determined by propagation require- 
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ments. This antenna must also comply with other requirements 
which insure complete transmission of electromagnetic energy 
from the generator to free space, when transmitting, or from 
free space to the receiver, when receiving. The problem is to 
find an adequate configuration and material for this device 
in order to verify the previously mentioned requirements. In 
general there is no routine method for solving this problem. 
Since the early days of electromagnetic waves attention has 
been directed to applying metallic conductors. Recently, loss- 
less dielectric materials were introduced [1] and lately, semi- 
conductors [2, 3} and superconductors [4] are being success- 
fully applied in design of active and passive antennas. 

Disregarding the problem of material choice, it is very 
difficult to find by analysis a metallic antenna configuration 
which complies with given requirements of impedance and 
directivity. Development in this area is a matter of experience 
and convenience. It is striking to mention that when the com- 
puter was invented — with its great storage capacity — some 
workers in this field thought that it could replace the require- 
ment of experimental work as it is able to give the needed 
information, just by solving Maxwell's equations with boundary 
conditions on the antenna surface [S]. In this statement though, 
no attention was given to how to choose this surface. The 
antenna problem however, is such that neither the current 
nor the surface configuration are known. Classical antenna 
design techniques, or antenna theory, is based on choice of 
the configuration and then finding the current distribution. 
the impedance, the radiation pattern and the frequency cha- 
racteristics. The difficulty of choice of the configuration to suit 
required radiation patterns has stimulated modern develop- 
ments of antenna arrays [6], where a lot of flexibility exists in 
arranging and feeding different elements. 

The object of this paper is to throw some light on usual 
antenna design techniques and a trial to find an adequate way 
which can be followed in the design considerations. As the 
subject is too comprehensive to be adequately eondensed into 
a small article, few points are carefully selected to give clear 
representation, without loss of generality, applicable to most 
antenna types. 


2. General considerations 


Developments of antenna theory are always mentioned in 
connection with solutions of Maxwell's equations for different 
known antenna configurations, and a classical discussion of 
antennas is usually based on some simplified assumptions 
which do not fit very exactly into a rigorous theory. This began 
in the early days of electromagnetic wave theory with the 
discussion of Hertz's experiments. Hertz is the founder of 
antenna theory; after few years from the publication of Max- 
well’s treatise on electromagnetic theory he has solved the first. 
the simplest and the most basic antenna problem of a small 
capacitor antenna (oscillating doublet or Hertzian dipole). 
After the year 1900 [7, 8] a number of very important theoretical 
papers were published on spherical, ellipsoidal, cylindrical and 
biconical antennas. These gave better understanding of the 
phenomenon. For time dependence of el” the solution of Max- 
well's equations (with the boundary conditions on the surface 
of the antenna) is found possible only for certain values of 
(the so-called proper values) and these are derived as complex 
quantities. From these complex quantities the O-factor of the 
antenna can be defined by analogy to normal resonance circuits. 
fw = @, + jw, from the free vibration of the antenna 0, will 
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be @,/2w,. It has been learned also that the Q-factor of spherical 
antennas for different modes are very small. some of which are 
even zero, which corresponds to damping without oscillations. 
This certainly showed that spherical antennas are typical ex- 
amples of wide-band antennas. 


3. Linear antennas 


The current along a thin ellipsoidal antenna and long prolate 
spheroids shows practically sinusoidal distribution when the 
length of the antenna is half a wavelength. which corresponds 
to resonance. The distribution deviates from this sinusoidal 
form away from resonance. The sinusoidal form is the base 
for what is called the conventional theory of antennas. Linear 
antennas of half a wavelength are assumed, with a good degree 
of approximation, to carry currents having this form. 

More accurate determination of the form of this distribu- 
tion is required for longer antennas. Trials to achieve this re- 
quirement began by Hallén for getting an approximate solu- 
tion for the z-component of the vector potential, A, due 
to currents flowing in the longitudinal direction of a linear 
antenna parallel to the z-axis. The differential equation for 
this vector potential is given bij 


ö Ĳ' 
— +k |A (2) =-— jj V, à (z) 
dz” WO 
where k = 2r/d, À being the wavelength of electromagnetic 


radiation; Ò(z) is the Dirac delta-function and V, is the applied 
voltage at the central gap of the linear antenna, whose length to 
diameter ratio is assumed to be large. The solution of this 
differential equation yields the following integral equation for 
the current distribution [9] 


he 


eikr / me ( 
Í 1_(z’) E dz’ — 4 jAri(u le) +4 C cos kz + 7 V, sin kiz 


h | 


where h is half the antenna height, r = (zz)? + a?}*, ais the 
radius of the conductor, ka <1, a hand C is a constant 
which can be determined from the boundary conditions. 
King [9] and his colleagues proposed different methods for 
solving this integral equation with some approximation. They 
could develop lately various practical solutions for this integral 
equation. Their three-term method gave good approximation 
to practical measurements. They could also extend their five- 
term method, which is concerned with two parallel antennas to 
include any number of parallel antennas. A nice survey of these 
efforts is given by King, Wu and Chen in a recent publication [10]. 

The efforts directed to the solution of the previously men- 
toned integral equation in nearly 80 years did not contribute 
as much to the design techniques of linear antennas. as the 
contribution of the conventional theory of antennas which is 
based on the concept of the sinusoidal form of the current 
distribution. This simple concept and others which are intro- 
duced into the antenna theory early after 1900 could lead to 
a lot of developments in antenna design techniques. One of 
these developments, and the most important one, is the Yagi- 
Uda antenna, which has been in use since 1926. Some practical 
and theoretical design techniques of this antenna [11] appeared 
in 1954. More accurate calculations for the current distribution 
on this antenna have not been available till the last few years when 
many workers began to formulate the problem systernatically 
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by the availability of the previously mentioned expressions for 
multielement systems [12]. 

The bandwidth of a linear antenna is determined by its 
thickness. This phenomenon has been learnt from the charac- 
teristics of the prolate spheroid. It has been proved that the 
thinner the prolate spheroid in question, the higher the equi- 
valent Q and the smaller the frequency bandwidth. Very thick 
antennas approach, in their configuration, the previously men- 
tioned sphere, which has been proved to be a typical example 
of a broadband antenna. 


4. Frequency-independent antennas 


These are known since 1957 when Rumsey broke the antenna 
bandwidth barrier with his angle concept. This concept was 
already known to Rumsey earlier than 1957, but he delayed 
its publication until Dyson demonstrated in September 1957, 
in the antenna laboratory of Illinois University, that the log- 
spiral is indeed frequency-independent. The representation 
of the field generated by a frequency-independent antenna is 
easy if Maxwell's equations are written in a form for which 
the differentiations are performed with respect to the normalized 
dimensions to the wavelength. When these differentiations are 
represented by rot’, Maxwell’s equations reduce to [13]: 


rot E oc Handrot H ax E, 


where the constant of proportionality is independent of the 
frequency. This means that if the tangential field on the surface 
of the antenna is fixed, and the dimensions of the antenna 
scale with the wavelength, the antenna will be frequency- 
independent. The next problem now is how to find the antenna 
which satisfies this specification and which suits the required 
application. At this point the ability of the human brain to 
understand the physical phenomenon and find out the adequate 
configuration of the antenna comes into action. In 1957 Du- 
Hamel and Isbell invented the log-periodic antenna which has 
found more wide applications than the original antenna de- 
scribed by Rumsey. It is interesting to mention also that some 
designs of the log-periodic antenna which were first proposed 
by Turner did not work (these designs were not published but 
remain in the history of development of these antennas). In 
1959 Isbell’s first log-periodic antenna could work. Lo [14] 
grasped the concept of an active region which moves with 
frequency and saw the necessity of the twisted feeder or reversal 
of alternate dipoles. This concept and the universal relation 
which relates the antenna impedance to the impedance of its 
complementary structure (4Z,Z, — Z°) are the most power- 
ful tools in this field. The analysis of these antennas (some 
workers in this field call it ‘theoretical work’; to my idea the 
theoretical work includes also the stage of finding out the 
shape of the antenna structure as well as its material) lagged 
in time, until experimental results were available for inter- 
pretation, because a rigorous solution of Maxwells equations 
is prohibitively difficult for their complicated structure. Mittra 
and his colleagues published some papers in the last few years 
on these analyses [15, 16]. 


5. Microwave and superdirective antennas 
In the field of microwave antennas the developing procedure 
did not differ much from that of the previously mentioned 


types of antennas except that the idea of using metalhiec reflectors 
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and lenses as antennas is taken from the known optical devices, 
The most useful concept which is of wide application in this 
field is the fourier transform relation between the radiation 
pattern of an antenna and its aperture field distribution [17]. 
The author succeeded to find mathematical concepts which 
were helpful in simplifying the design techniques of parabolic 
and other modified types of reflector antennas. The most 
simple concept is what is called the edge relation for the illumina- 
tion function at the reflector edge [18, p. 46). 

In the field of superdirective antennas many efforts are in 
process since the publication of Bouwkamp and de Bruijn [19] 
in which they have shown that there is no limit for increasing 
the directivity of the antenna by choice of an idealized current 
distribution. There are practical limitations for this phenomenon 
because of the large amplitudes and high rates of phase varia- 
tion of these currents on the antenna structure [20]. This un- 
favourable condition arises not only when it is attempted to 
increase the directive gain excessively, but also when it is re- 
quired to closely approach a desired nonrealizable directive 
pattern or to realize unjustified decrease in the sidelobe level 
(21, p. 1788). An example of the phenomenon of superdirec- 
tivity which is associated with excessive decrease in the side 
lobe level is given in [18, p. 92]. 


6. Conclusions 


From the previous discussion it is easy to see that the antenna 
design techniques present the following three stages: 

|. The stage of finding an adequate configuration of the an- 
tenna to satisfy the radiation requirements and to suit the 
field of application. In this stage idealized situations can be 
studied, from which simplified assumptions and conclusions 
may be obtained. This stage is characterized by fascinating 
simplification of the electromagnetic theory in combination 
with observation which give physical insight into the problem. 
Successful results in this stage depend on our stored know- 
ledge as to the behaviour of fields and currents. This can 
help in correlating the simple facts of electromagnetic phe- 
nomena and constructing a pattern for the current distribution 
and its supporting structure in order to arrive at the adequate 
configuration. It is interesting to mention here the words of 
Rumsey in the introduction of his book on frequency-inde- 
pendent antennas [13]: ‘Research in this area has uncovered 
many striking phenomena which remain unexplained, and has 
generated solutions of Maxwells equations which have yet to 
be observed’. 


2. In most of the cases the proposed construction in the previous 
stage may be complicated or it may not work. Therefore it is 
preferred preliminarily to test the possible practical applica- 
bility of these antennas before proceeding to the stage of anal- 
ysis. This helps, firstly to gain more experience for inter- 
pretation and guidance of the assumptions for mathematical 
solutions of the complicated configurations, and secondly to 
avoid the risk of time waste if, for some reason or another 
which could not be seen in the first stage, the originally pro- 
posed constructions are not practically realizable. This stage, 
and the previous stage, are especially essential for antennas in 
which unconventional matertals, such as semiconductors, may 
be used. 

In addition to this conclusion it is interesting to mention 
these words of Levi [22]: ‘Measurement is fundamental to all 
science, and the formulation of concepts used in such measure- 
ment is itself an important part of the science’. 
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3. Experience from the previous two stages is essential for 
performing useful mathematical analysis. The object of such 
analysis is to help the designer in evolving his original physical 
picture and bringing order to the multiplicity of facts obtained 
in the first two stages. When, sometimes, a quantitative solution 
of the problem is hard to find, those physical pictures often give 
qualitative answers. Maxwell’s equations themselves express a 
physical point of view. At this point it is interesting to mention 
the definition of ‘practical solution’ as given by King [9, Chapt. 
9): ‘Solutions are called practical if they combine reasonable 
quantitative accuracy with formal simplicity in a physically 
meaningful representation. Such solutions can be used. on the 
one hand, to solve engineering problems and, on the other 
hand, to illuminate the underlying physical principles in their 
relation to the actual situation’. 

As the simplicity of the mathematical formulation of an- 
tenna problems contrasts sharply with the complexity of its 
solution, the stage of developing a useful practical solution is 
the longest and difficult stage, especially when the designer is 


guided by others’ physical concepts. Some designers, however, ° 


use concepts which are inconsistent with Maxwells equations 
and yet with this aid some approximate results have been ob- 
tained. Such a concept is that the power losses due to radiation 
produce the same effect on transmission of waves along con- 
ducting wires as losses due to dissipation. We hope that the 
use of modern computers with their modern useful applications 
will give a hand to help shortening this stage of antenna design 
techniques. 

Finally, it is clear to see that the theory of antennas with its 
three closely related stages is a typical example of the close 
relation between science and technology, the relation described 
by Bondi [23] as that of ‘the chicken and the egg’: you cannot 
have the one without the other. This is actually the physical 
nature of science. 
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1. Introduction 


In considerations regarding the formation of images through 
the illumination of a hologram, authors always use, from the 
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Asymptotic theory of holographic image formation 
by dr. ir. H. J. Frankena, Detft University of Technology, Department of Electrical Engineering 
Summary: Using the exact formulas of scalar diffraction theory, the field behind an illuminated hologram is 


investigated. The results are expanded asymptotically as k — oo, leadin g to expressions which describe well-known 
phenomena that usually are deduced from approximate expressions valid only for restricted classes of configura- 


smallness of the optical wavelength as compared to the spatial 
distances in the configuration and from the smallness of the 
object, the hologram and the images as compared to the dis- 
tances between them. The present account intends to deliver 
a more exact theory on the same subject, thereby only em- 
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phasizing the smallness of the wavelength with respect to the 
various distances in the configuration. Our aim is to show 
the application of certain methods rather than to give the most 
general proofs. Therefore, we often shall restrict the account 
to simple cases. 

The electromagnetic field will be represented by the complex 
scalar function U (x,y,z;0) (see, e.g., [1]), where we assume the 
field to vary harmonically in time with the factor exp( —imt), 
this factor is suppressed in all subsequent formulas. The func- 
tion U is a solution of Helmholtz’s equation 


VEU + kEU = 0, (1.1) 
where k = w/c is the wave number and c« = velocity of propaga- 
tion of the optical waves. : 


We follow the well-known procedure in which we firstly 
calculate the intensity | U|* at the place of the hologram during 
its formation and then, assuming the transparency of the holo- 
gram after development to be proportional to this intensity, 
calculate the field created by a wave after its transmission 
through the hologram during the reconstruction phase. 

The field U during the formation phase consists of the super- 
position of the field U, coming from the illuminated object 
and of the field U, usually called the reference wave. From 
a known object, illuminated by a field U of angular frequency 
_w, and placed in a homogeneous and isotropic medium we have 
the contribution [2]: 


U» (x, y,2, w,) a || je (ns Von Eb w‚)jn Ld 
ola jee 1 


exp(ik, R‚) 
drnR 


al 


_U ons Van: Sob’ wk, 15 Rs) on,jnf d S ob 


(1.2) 


where d/Ôn denotes the derivative at the surface of the object 
in the direction of the inward normal and k, is the wave number 
in this situation. Here, R, stands for 


Rav Dn (x Ki) je (y — Vo 5 E (z Ee, Zoo) ; ì s Ö (1 3) 


Choosing the medium to be lossless we know that k, is real. 
For the reference wave we insert a plane wave of the same time 
dependence as U. travelling in a direction given by the phase 
vector fl. Its field U, is then 


Urer (x, NY, Z. w,) zen A, Exp B, 7) (1.4) 
where 
Pa xl, + yÂ, + zï (1.5) 


and A, is a constant. We imagine the hologram to be situated 
in the plane z = 0 and the fields U, and U, are assumed to 
be incident from the half-space z < 0. The total intensity of 
light / incident upon the film is then 


1x, y) =| Up + Urerl Ee or 
=| 0 0) + | Ur Go. OT DI + 
+ Ui ey, 0 5 0) Uber P,O 50) + 
+ UE (FO 500) ier Cy. 0 5 0) (1.6) 


We suppose, in accordance with statements from specialists 
in this field, that the transparency t of the developed hologram 
is directly proportional to /. If, during reconstruction, a wave 
coming from z < 0 strikes the hologram, the field transmitted 
through the hologram will be 
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U (x,y, 0” U) = (x,y) Urs (x, Y, 0: 0) 
(x,y on the hologram). (1.7) 


where U, is the wave function describing that wave. We 
assume the angular frequency during reconstruction to be @,. 
Taking for U again a plane wave, we have 


Weis (x, Vz, @) IK A; Exp ip, 7) (1.8) 


where A, is a constant. With (1.7), the function U(x, y, 0” ;w,) 
is, apart from a constant factor, equal to the product of the 
expressions at the right-hand sides of (1.6) and (1.8). At points 
in the half-space z > 0, the field U can be obtained by application 
of a representation theorem. We let the hologram occupy the 
rectangle —aSxSa, —-b<ysb, z=0; outside this rectangle 
we take the factor f in (1.7) to be zero. Let the wave number in 
this situation be k,; then, from the representation theorem for 
a half-space [3] 


aû rij 


dl io Vim’ Zim’ w.) mice a | | U (x, y, 0; w) X 
| An 
exp (ik, R 
0 JP Kie) da (1.9) 
ÓZin Kn 

we arrive with (1.5)... (1.7) at 
a h 

U is Vin ein” @,) Sik [ | Ue (x, y, 0; w) Xx 
_ —h 


exp (ik, R‚) 
f ee U f a | Dsl E Us Wi E 


bd eN 
Ozi Kn 
+US U ‚| dxdy (z. > 0) (1.10) 
In these formulas, R‚, denotes 
Rin = [in 2)? + Oi VIE + Zi 1° 20 (1.11) 


c is a constant factor while the arguments of U, and U; in 
the expression between square brackets have been given in (1.6). 


2. Discussion, of the various contributions to the wave function 
during reconstruction 


In order to discuss the result (1.10), we consider the contributions 
of each term between the square brackets separately. 


a. Inserting (1.8), we obtain from the first term a contribution 
of the general form 


a h 
Ö 
U, = F | [ L&P exp (if, ,X+ IB, V+ 
im —a —b 
+ik, Ri) Ri, dxdy (2.1) 


where in this case f, (x, y) contains information concerning the 
amplitude of the field coming from the object but not concerning 
its phase. In many calculations, authors do neglect this term 
in comparison with the next term in (1.10). This procedure 
seems to be incorrect, as most experimentors try to obtain almost 
equal values of |U, |* and | Ul“ at the hologram plane during 
the formation. | 

In view of the fact that k, has a large value for optical fre- 
quencies, we can expand the integral in (2.1) asymptotically 
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as k,—c0. Thereby, we assume | Ul* to vary slowly as com- 
pared with the exponential functions in (2.1) (both B, ‚ and 
fB, „contain a factor k,). We apply the method of stationary 
phase [4]; for simplicity, we restrict the considerations to the 
first approximation in which only the contribution of the 
stationary point is taken into account. Introducing the direction 


COSINES Y, … 7» „and 7, ‚ through 


B, sz Kala B, ,= Kaas 


8 


nn kmi K B 2 
Ln De “Nas 420 


we obtain from (2.1) by subsequent application of the method 
of stationary phase on the integrations with respect to vand to x: 
U, zr F, irm» Viens Zim) Xx (2.3) 
x EXp ik, (ss Lim of Van Fin ij 12e Em) (Zin =0, kh, 00) 


where F, is an algebraic function of the coordinates En 
Zin in which a.o, the function f, from (2.1), calculated in the 
stationary point x = x,, Vasil Ya zl V S Vins —V2,y| Zina. 0 | 
appears as a factor. We do not consider the special case Pz. =0 
(striking incidence on the hologram). which has no practical use. 
From this, we obtain a non-zero contribution for ANY Xine Vins 
Zim in the relevant half-space; it forms nowhere an image and 
acts as a spurious illumination, spoiling possible \mages created 
by the further contributions to U in (1.10). We remark that 
the asymptotic expansion can be extended to one giving con- 
tributions of higher-order terms too, such as the contributions 
from the end points of the intervals of integration. However, 
such an extension will not lead to a solution that yields a better 
understanding of the phenomenon described by U. 


b. Considering the second term between square brackets in 
(1.10) we remark that, from (1.4), U erl° is only a constant 
while U (x, yv, z; w,) is a solution of Helmholtz's equation 
with the same angular frequency wW, as occurs in the factor 
between braces in (1.10). Then, on account of the uniqueness 
of the contribution U, as given by 


exp(ik, R 


im) dxdy 
R 


a h 
U,=c| Al} Í Í U. (x, poe) 
“im im 
-a —b 


(2.4) 


for a known distribution Uice (X% y, 0; ww), the field EE 18, 
apart from a constant factor, the plane wave U from (1.8) 
which propagates through a rectangular aperture into the half- 
space z, > 0. 

Hence, this second term describes a contribution which 
consists of the reconstructing plane wave which travels through 
the hologram into the half-space z,, > 0, thereby undergoing 
a change in amplitude. 


c. Inserting (1.2), (1.4) and (1.8) into the third term of (1.10), 
we arrive at the contribution U,: 


Be 
Özin abject K: 


h 


Jz (ons Van? Zob' XF, 0) Ri, Rt x 
h 


Xx EXP (ik, Ro, + ik, ya „Xx + ik, pop + ik, Rumi Vie X— 


— ik, 7,,) dxdy | dS, (2.5) 


| 
| 


+ 
where 7, and Y1,, come from f#, in a way similar to (2.2). 
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Here, f, is the function between braces in (1.2) multiplied by 
a constant. For simplicity, we restrict the calculations to the 
case k, = k, = k and, again, try to obtain from (2.5) an asymp- 
totic expansion for U, as k — oo by means of the method of 
stationary phase. Furthermore, we put Vi,x S Vz, x Ad, = 
72, ‚ Which means that U, and U, propagate in the same 
direction. 

Writing the exponential function at the right-hand side of 
(2.5) in the form exp { ikh; (x, v) }, we know that in the stationary 
point êh‚/ôx = O and £h‚/ôy = 0 simultaneously, or 
Ark) Ran = Cin AVR 


WV) Rop = Wim IN Rim 


In order to visualize the first relation, we consider Fig. l and 
remark that in the stationary point we have sin(a) = sin(a) 
or, in this case, x= «, (in Fig. 1, the points P , and P, have 
the co-ordinates x, yv, z, and Kim Vim Zims Fespectively). 
The set of points in the plane z = 0, satisfying this relation, is 
a curve of which in any case the point P (being the point at 
z = 0 of the straight line between P,, and P‚). is a member. 
Analogously, the second relation (2.6) yields a set of points 
containing P. Hence, we may conclude that P is the stationary 
point; its co-ordinates are 


(2.6) 


onZim — Kin Zord Zien Zg) Wes Zim —YimZob)l Zim — Zoi), 0 
(recall that z, < 0). Further we remark that for k, = k, we 
have 
hz (x,y) = Ro + Rim (2.7) 
which in the stationary point P is equal to the distance Rasta 
between the points P, and P… 
After application of the method of stationary phase at the 
integrals with respect to x and y, we obtain from (2.5) 
. explik R) 

Ï pe 5 me an : E ali irm 

b 3 Í F, (A op: } ab" ST uh* Kien’ } im” Zin) De Wan sit ds 


un je akai rm 


ah 


(hk — ) (2.8) 


in which F, is, apart from a constant factor. equal to the product 
of the expression between braces in (1.2) (which itself. provided 


Fig. 1. The geometry from which the stationary point P is deduced. 
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that R‚, # 0 (object outside hologram plane), can be reduced 
for k — «o) and an algebraic function of the co-ordinates of 
both P, and P,…. This function, which is related to the second 
derivatives of h, in the stationary point, will be almost constant if 
2, varies only slightly over the object and if |z,, is large as 
compared with |; — Xe, | and [Vip op |- In that case, the right- 
hand side of (2.8) is proportional to the right-hand side of (1.2) 
with x, y, z replaced by Xin» Vins Zij: This means that then the 
field U, for zi, > O is, apart from a constant factor, the analytic 
continuation in z,, > 0 of the field which was created by the 
object in its original position. It means the occurrence of a 
‘virtual image’ at the same place as where the object was situated 
during the formation phase. If, however, |z| is not large 
enough, the algebraic factor in F, will cause deformation of 
this image. 

The present considerations apply only to the case where the 
size and position of the hologram are such that the stationary 
point is inside the intervals —a<x <a, —b<y<b. In any other 
case, an asymptotic expansion will yield only the higher order 
effects due to the boundaries of the hologram and no virtual 
image can be seen. This is easily demonstrated in an experimental 
setup by choosing the position of an observer (Xin: Jims Zim in the 
formulas) improperly. 


d. The last term between square brackets in (1.10) yields the 
contribution 


„a È 


5 Û | | in rele 
Us= ek | Ĳ fu os Von: Ev XY, 0 ) Ro R, x 
OZ object | | 


x expl—ik, R, + ik, Yan + ikea ,V + ik, Rim + TE 


+ ik, 1,49) aud) dS, (2.9) 


with the same notation as used in (2.5). The function f, differs, 
apart from constants, from f, in the complex conjugation of the 
expression between braces at the right-hand side of (1.2). Again, 
we put k, = k, = k and we write the exponential factor in the 
integrand as exp {ikh, (x, v)}. 

The considerations are now restricted to the simple case 
Pi S Pip S Tae = Y2,y = 0, Le. the case that both Us 
and U, are travelling along the z-axis. Then we have 


h‚ (k‚p) = Rs, + Kin (2.10) 
The stationary point follows from 
— Xn) Rop = in Xx) Rim (2.11) 


Wy Ken == HSN Rim 


It is now simpler to proceed algebraically than to follow the 
geometrical observations of the preceding section. 
From (2.11) we deduce (as z,, < 0) 
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(AX) Ain) SWV) Ii) S Ron Rim 2 ov/ Zim (2-12) 
From (2.10) we can write 


hem ER? + Ra? ZR Rim) 13) 


in which, in the stationary point, with (2.12) 
Rin Ran = Kim CHA) = | 

= (XK) Kim) + or) Wima) — Zot Zim (2-14) 
Insertion of (1.3) (with z = 0), (L.1L) and (2.14) in (2.13) leads to 
he= Fo Xin)® + Via) + (Zon +Zim) 5 ERovim (2-15) 


in the stationary point, where Rij, is the distance from the 
observation point Xi, Vins Zim to the point Xp, Yab: — Zop: The 
sign in (2.15) has to be chosen in accordance with (2.10), in 
which R, and Ri, are both positive. Hence, in the domain 
Ri, Z Ri, we arrive at an asymptotic expansion for the inte- 
grations with respect to x and y in (2.9) of the form exp (ik Rovi) 
times an algebraic factor which is practically constant when 
|z| varies only slightly over the object, while |z| > | Xin — ob! 
and [ziel > Win Yool. For Rim < Roos this expansion has 
exp (—ik Rij) as the exponential factor. Now, we obtain for 
U. 


ï ns Û F lei bl 
L Et | | Fi Cop: Vobs Zobs Fims Vime Zim) à 


* © geet 


exp Gik R… 
xp CI B min) d8, 


abiim 


(k — 2} (2.16) 


XE 


For R‚, < Raj this means waves converging to a domain found 
by mirroring the object with respect to the plane z = 0; for 
R, 2 Ro, these waves are diverging from that domain. This 
means, physically, that a ‘real image’ is formed there. All remarks 
in section 2c. concerning the deformation of the image by the 
above-mentioned algebraic factor and concerning the non- 
existence of this image if the size and the place of the hologram 
are ill-chosen, apply here too. 
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Een gegeneraliseerde Schlömilch-reeks 


door ir. H. W. F. van 't Groenewout, rotterdom-Overschie 
Dedicated to Professor J. P. Schouten on the occasion of his retirement 


Summary: A generalized Schlämilch series. 


An itemized survey is given of several infinite Bessel series which arc closely related to the well-known Schlömilch 
series. The series studied were derived from an investigation of the quantization noise that goes with a pulse-code- 


modulation system. 


1. Inleiding 


Toen kort na de oorlog, mede onder invloed van het baan- 
brekende werk van Shannon op het gebied van de informatie- 
theorie, belangstelling ontstond voor nieuwe modulatiesystemen, 
gaf Schouten een zeer oorspronkelijke methode aan om de ver- 
vorming te berekenen die bij pulscode-modulatie (PCM) in de 
vorm van quantiseringsruis optreedt. Daarbij werd gebruik ge- 
maakt van de laplacetransformatie om de sprongen die ontstaan 
als het signaal achtereenvolgens de equidistante drempels over- 
schrijdt, mathematisch weer te geven. De resultaten van deze 
studie werden gepubliceerd en naast een inzicht in het juiste 
gedrag van de modulatieprodukten werd een bevestiging ge- 
vonden van de grootte van de signaal-ruisverhouding van een 
dergelijk modulatiesysteem zoals die door anderen langs geheel 
verschillende weg was gevonden [1, 2]. 

Bij deze studie kwam het, kort samengevat, hier op neer dat 
de gezochte quantiseringsruis kon worden uitgedrukt als de 
fourierontwikkeling van de zaagtandfunctie waarbij als abscis 
niet de tijd zelve doch het ingangssignaal als functie van de tijd 
gekozen wordt. Dit leidt voor bijv. sinusvormige ingangssig- 
nalen tot uitdrukkingen van het type 


Ens et Ï s ï | 
N sin (Zr A sin (w IN 1 
IT 

zr 
Een nadere uitwerking van deze reeksontwikkeling leert dat de 
amplitudes van de gezochte modulatiecomponenten gegeven 
worden door reeksen van het type 


id ad (2rrz) (2) 


een oneindig voortlopende reeks dus van Besselfuncties van 
de eerste soort waarbij het argument van iedere term evenredig 
is met de rang van de term. Reeksen van dit type zijn eerder 
onderzocht door Schlömilch en anderen, doch voor het be- 
trokken onderzoek bleek het noodzakelijk te zijn een verder- 
gaande studie te maken van reeksen van een meer algemeen 
karakter en wel 


at 


L =| ad, (@zrz) pl2zry) (3) 


waarin p(2rry) een goniometrische functie voorstelt. 

Gezien de bijzonder snelle en elegante wijze waarop men met 
de laplacetransformatie het doel had bereikt, lag het voor de 
hand om ook voor de sommering van deze bijzondere reeksen 
van dezelfde transformatiemethode gebruik te maken. De 
resultaten van de berekeningen langs deze weg en wel van 
reeksen van het type: 


Jas, Crrz) 


(4) 


Hee 
de | 


(rerz) 
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werden in een intern rapport [3] vastgelegd en de toegepaste 
rekenmethode in een publikatie [1] in hoofdlijnen weergegeven. 


In het onderstaande zullen wij ons bezighouden met de sommatie 
van de meer algemene reeks 


Ke J, Fr (@rrz) 


Ë 
mij | 


cos (2rry), (7) 


(rerz) 


die behoort tot de categorie onder (3) genoemd en waarvan in 
wezen de reeksen (4), (5) en (6) als limietgevallen beschouwd 
kunnen worden voor respectievelijk y = O en y =— 5 terwijl de 
reeks (6) gevonden wordt door het verschil te nemen van (4) 
en (5). Bij deze berekening wordt om de mathematische aar- 
zelingen te overwinnen, die nu eenmaal inherent zijn aan een 
rechtstreekse toepassing van de laplacetransformatie bij de 
termsgewijze integratie van oneindig voortlopende reeksen 
waarvan de mate van convergentie a priori niet vaststaat, een 
meer klassieke analytische weg gevolgd. 


2. Berekening voor n — 0 


Als inleiding tot de berekening van de reeks (7) onderzoeken wij 
allereerst de reeks gegeven door 


AN 4 J, (Arerz) 


Fil 
mi = | 


cos (2rry) (8) 
(rerz) 


waarin J, de Besselfunctie van de eerste soort en orde v voorstelt 
en zowel v, y als z reëel verondersteld zijn. Voor z beperken we 
ons tot waarden groter dan nul, terwijl op grond van de periodi- 
citeit van de goniometrische functie zonder verlies van alge- 
meenheid y beperkt kan blijven tot het interval 0 < y < & 

Daartoe bepalen we allereerst de waarde van de uitdrukking 


NL J_ (2rrz) 


mais (rerz)' 


cos (arry) (9) 


Op grond nu van de integraalbetrekking van Poisson voor de 
Besselfunctie van de eerste soort en orde v: 


rÌ 


cos {x sin (d)}cos” (@) dep 
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A 


geldig voor v > — ; kan voor (9) geschreven worden: 


J (2rrz 
al Á did cos (Arry) =— wen 


1e en EET 
7 (rrz) 11 (» 5 7 (- ) 
2, 2 
x [: ln cos {2rrz sin (@)} cos (2rry) cos°° (@) dp (10) 


Na sommering van de eindige reeks onder het integraalteken 
en onder gebruikmaking van de bekende relatie voor de 11- 
functies: 


IH (xD) HI (y —I) 


=d IN cos! (p) sin°” '(#) dp, (11) 
H (x+y-l) 


ü 


geldig voor x > Oen y > 0, ofwel v > —5, gaat (10) over in: 


Es 1, Gn) cos (2zry) + a je 
TE Kezr 211») 211 pe H ze 
2) \ 2 
| 
| > sin{ORsDrlzsin 4) ze (5) as 
ö sin {mx (z sin (@) + y 


ER \ sin{2RyDe(zsin(g) IN or (9) a (12) 


0 sin {rr (z sin (@) — y) 


Wijziging van de integratievariabele in de eerste integraal van 
(12) door tf = zsin(d) + y te stellen en in de tweede integraal 
door de substitutie f = z sin(d)—y en overgang tot de limiet 
R — oo geeft als resultaat 


En sin ((2R +1) zr) - 
6 sin (st) 


br Sem * : ij liek 
+ | sn@RtDel 2 4 za (13) 
ä sin (rf) 


Afhankelijk nu van de waarden van (z + yv} en (z — y) kunnen 
drie gevallen onderscheiden worden: 


a. O<z<ysj.Inditgevalis0 <(z + y) < len —}<(z2—y) 
< 0, zodat door tekenwijziging in de tweede integraal voor de 
som van beide integralen geschreven kan worden: 
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RIE ANR 
(rj) (5) 
al \ 2 


Ie sin (2R+I) rr} {z° ds dt (14) 


s sin (st) 
Nu leert het onderzoek van de integraal 


le, td E 15) 


sin (sf) 


dat deze na substitutie van een nieuwe variabele 


Ë = (t—ytz)}/ly tzt) evenredig is met 


EN, 
a pr 


Í Ee (16) 
Ja (E+1) sin [r {y +(& — 1) z/(&+ DI] 


waarin de integrand = O(& "> als y # zen = O(& "7 d) 
als y = z. De integraaluitdrukking (16) convergeert dus absoluut 
alsv > — jeny# zenalsv > jen y = z, zodat op grond van 
het theorema van Riemann-Lebesque de integraaluitdrukking 
(14) = O(1). Daarmede vereenvoudigt de uitdrukking (13) zich 
tot het eindresultaat 
 J(Arrz) Í 
EE ATD 
(rrz) 2E) 


1 1 7 1 et 1 1 
als-; <vSjzen0<z<ysjenalsv>jen0<zsysi. 


=0 (17) 


r=l 


b. OSy<zstInditgevalis0 < (z + y) < len0 <(z—y) 
< }en beschouwing van de beide Dirichlet-integralen van de 
tweede soort in de uitdrukking (13) leert dat gezien het monotone 
verloop van de functie 

| 


Eelt 2 (18) 


voor R — oo de waarde van de integraal gelijk is aan 


Ì 
2 


| "-= 
ie tp eri lanl ‘fd: 
2 t=û 2 | if = Û 
waarmede (13) overgaat in het eindresultaat 


J (2rrz) | 
ee 608 (Zar y) ze 
(rerz) 211») 


En 


mp =| 


ENVER (19) 


alsv> —jen0Sy<zst. 

c. OS ysSis<z. Indit geval is(z + y) > len (z-y) > À. 
Stelnudat (z + yv) begrensd wordt door het interval k < (z4-y)< 
(k + 1), waarin k een geheel en positief getal (k = 1, 2, enz.) 
en (zy) begrensd wordt door het intervalm < (zy) < (m1), 
waarin m eveneens een geheel en positief getal (mm = 1, 2, enz.). 
Beschouwt men nu vervolgens de eerste Dirichletintegraal, dan 
ziet men dat de integratievariabele t alle gehele waarden van 
l tot en met k doorloopt. Splitst men nu vervolgens de integraal 
in de intervallen 0 — 1, 1 > 2, .…, (kl) > k, k — (k+1), 
dan levert iedere deelintegraal een bijdrage gelijk aan de halve 
som van de waarden van de functie (18) ter plaatse van de 
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ondergrens en de bovengrens, en het totaal wordt gegeven door 
de eindige reeks 


EN Kes 


bs 2 rs k 2 rl | 

la IP) Pt Der kg (20) 

Past men nu eenzelfde beschouwing toe op de tweede Dirichlet- 

integraal van (13) en sommeert men deze resultaten met die 

gevonden onder (20) dan luidt het gezochte eindresultaat : 
J,(@2rrz) î 


pT Cos (2rry) + — 
rel (zefz) 2) 


| dE 
tn | IN -Ek 


af) 


+ Pk ZA (ZP -mpf} d (21) 


als v > —jen k en m gehele positieve getallen, zó gekozen dat 
k < (zy) < (ktl)enm < (zy) < (m+l). 


3. Discussie 


De bovenstaande resultaten, samengevat in de vergelijkingen 
(17), (19) en (21) leiden tot een groot aantal interessante resul- 


taten voor zekere waarden van z en y waarvan enkele in de 

literatuur bekend zijn. Zo levert (17) voor v = Ì de reeks [7]: 
ES ‚ d, (&rrz) | 

LN + = (22) 

û (erz)” 211) 

als —<vSjen0 <z<jenalsy > ben0<z 5 Hieruit 

vloeit als bijzonder geval voort de reeks: 
ED | | 

Db (IJ, (arrz) 4 =0 (23) 
r=l p) 


geldig voor 0 < z < & [4]. 
Voor waarden van z, groter dan ; geldt, indien y = &, de 
betrekking k = (m1), waarmede (21) overgaat in: 


J (2rrz) | l 


ap F F 
UO Etn 
= v— | 


(rerz)’ 21I(v) 


: In I\? rg | 
EN ml De dj (24) 


voorz > (en wel zodanig dat (n —5) < 2 < (n + 1), een resultaat 
dat eerder gevonden werd door Nielsen [6]. 
Stelt men daarentegen in (19) y = 0, dan ontstaat de reeks: 


» J,‚(@2rrz) | l 


Vi ner HMH 
id id ens) 
2) \ 2 


geldig voor v > —jen0 <zs 5, een resultaat dat eveneens 
van Nielsen afkomstig is [6]. 
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4. Berekening voor » > 0 


Tot slot kan nu de gevraagde reeks (7) uit de berekende reeks (8) 
afgeleid worden door op de linker- en rechterleden van de eind- 
resultaten (17), (19) en (21) n malen de operator 


ld/_d 
ken see) (26) 


toe te passen en vervolgens de macht van de noemer die door het 
herhaalde differentiëren tot v—2n gedaald is, weer gelijk aan v 
te stellen. Men verkrijgt dan uit (17) de reeks 


2 
—amr(d ld dan (2rry) =— (27) 
(rerz)' 


mm si | 


met geheel en > 0, geldigals — <vSten0sz<yst 


enalsv > jen0 <25 ys}, De reeks (19) levert vervolgens on 


(2rxrz) —)y 


„oe Jans CON dr pn 


eK 
(rra) gent, en (antr )A(—) 
| 2 2 


NAT 
pe \ an | (28) 


met n geheel en > 0, geldig als v > — jen0Sy<zsi 
terwijl tot slot de reeks (21) als resultaat oplevert: 


„J _)" 
mn el bede 


PELT Ona Wl = 
2 2 


| 
Andr — 


(k yy | à È (29) 


Eu | Í ( d À fi 1 (rm + y) ak 
mt = | (em + py)" dz") \ z° 


met » geheel en > 0, geldig als v > — jen k en m gehele po- 
sitieve getallen zo gekozen, dat k < (z + v) < (k + 1) en 
ms (zy) < (m1). De functies in het rechterlid, hier weerge- 
geven in operatorvorm, kunnen herleid worden tot vele bekende 
functies, zoals die van Gegenbauer in het meest algemene geval, 
polynomen van Legendre voor v = 5, en polynomen van Tcheby- 
cheff voor v = 1, terwijl voorts ook asymptotische benaderingen 
voor grote waarden van z gevonden kunnen worden. Het zou 
echter te ver voeren om deze mogelijkheden hier in detail te 
vermelden. 


k | EN: in Í 
4 mf „1 — ; 
‚ | (k yy" beed l z° | 


PF ] 


5. Slotbeschouwing 


Het eindresultaat, vervat in (29) kan ook — zij het in iets andere 
vorm — gegeven worden door niet met behulp van (26) naar de 
variabele z te differentiëren, doch door directe differentiatie 
van (17), (19) en (21) naar yv. Men verkrijgt dan bijv. uit (21) 
als resultaat: 
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A 


J,. (Arrz) 


N08 Arre 
€ 5 
(—) 
peren Dal (Sen 
| 2 2 dy. 
| Ï 
ok 3 E] de m2 ES kde 
+) le -(k-)| od Nen a | (30) 


onder dezelfde condities als bij (29) vermeld. 

Differentieert men niet 2n-maal, doch (2n— 1 )-maal (n= |, 2, 
enz.) dan ontstaan reeksen met p (2nry) = sin (2rry), waarvan 
door uitwerking van (17) als voorbeeld genoemd moge worden: 
Ss Jer (2rxrz) 


merel 


sin (2rery) — 0 (31) 
Grz) 
met n geheelen > 0,als —4<vsjen0Sz<ysSjenals 


259.51 


v=den0 Sz5SysS 5 Het is duidelijk dat ook hier weer een 
terugvoering naar vele bekende functies mogelijk is. 
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The surface line source problem in elastodynamics __n 


by prof. dr. ir. A. T. de Hoop, Delft University of Technology, 


Department of Electrical Engineering 


Summary: At the plane boundary of a semi-infinite, homogenebus, isotropic, perfectly elastic solid a pulsed line 
source of normal pressure generates a two-dimensional elastic wave motion in the solid. The displacement vector of 


this elastic wave motion is calculated at any point inside the solid. 


1. Introduction 


The surface line source problem in elastodynamics involves the 
calculation of the displacement vector of the elastic wave motion 
in a semi-infinite, elastic solid, generated by a pulsed line source 
of force applied at its plane boundary. Lamb [1], who has been 
the first to investigate this problem, calculated the displacement 


vector at the boundary. In the present paper the displacement. 


vector at any point inside the solid is determined by employing 
the author's modification (de Hoop [2)) of Cagniard’s technique 
for solving seismic pulse problems. 


2. Formulation of the problem 


The elastic waves under consideration are small-amplitude 
disturbances travelling in a semi-infinite, homogeneous, isotro- 
pic, perfectly elastic solid. The physical properties of the solid 
are characterized by its mass density p and its Lamé constants 
Â and u. A right-handed Cartesian co-ordinate system x, yv, z 
is introduced such that the elastic medium occupies the half- 
space —o<x<o0, —0<y<o,0<z<o0. A point inside 
the solid or at its boundary z = Ois located by either its Cartestan 
co-ordinates or its cylindrical co-ordinates defined through 


r cos((l) (2.1) 


Il 


x= r sin (4); y=}; z 


with0 Sr < oo, -—n/2 < 0 S n/2. The time co-ordinate is 
denoted by t. 

Along the line x=0, —co<y<oo, z=0 a pulsed force is 
applied to the free surface of the solid. It starts to act at the 


ELEKTRONICA EN TELECOMMUNICATIE 1 / 30 JANUARI 1970 


instant £ =0; prior to this instant the medium is assumed to 
be at rest (Fig. 1). 

Let d = Ï(x, yv, z, t) be the displacement vector of the elastic 
wave. At any interior point of the solid the stress tensor, 
with components Tt tT... T‚‚s iS related to the displace- 
ment vector through the (linearized) constitutive relations 


tr. = Adivù + 2u(ôu,/ôx), etc. (2.2) 
Ty = Te = HOU, /Ôy + Ou, /Ôx), etc. (2.3) 


From the (linearized) equation of motion and the constitutive 
relations it follows that, in the absence of body forces, the 
displacement vector satisfies the elastodynamic wave equation 


pp? grad div # — ve? rot rot ú — O?u/ôr® = Ô (2.4) 
mn which 
vp = [A + 24)/p)' (2.5) 


is the velocity of propagation of compressional or P-waves 
(for which rot # = Ó) and 


vs = (u/p)' (2.6) 
is the velocity of propagation of shear or S-waves (for which 
div u = 0). 

Now we investigate the case where the amplitude of the 
applied force is independent of y; then the boundary conditions 
as 210 are independent of y. Since the elastodynamic wave 
equation admits y-independent solutions the generated elastic 
wave motion is independent of y. Hence, #° = ix, z, f). The 
y-independent solutions of (2.4) separate into two classes: in 
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=Lrr = Eitdin) 
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“elastic solid 
Â‚p 


ür 


rd 


Fig. 1. The configuration to be investigated. 


one of them (has only a y-component, in the other 1 has only 
an x- and a z-component. Let the plane of the boundary be 
‘horizontal’, then the former displacement vector corresponds 
to an SH-wave (horizontally polarized shear wave) while the 
latter displacement vector corresponds to a superposition of a 
P-wave and an SV-wave (vertically polarized shear wave) 
(Bullen [3)). The calculation of the SH-wave is a simple problem 
in scalar wave propagation and is not considered here. The 
calculation of the combined P- and SV-wave motion is the 
subject of our further investigation. 

Let &” be the compressional part and #*Y the shear part 
of the relevant displacement vector, then ú = ùP + usw, 
From (2.4) it follows that @®’ satisfies the two- dimensional 
wave equation 


(O2/0x? + 82/02? — vo 20E = (2D 
together with the auxiliary relation 7 i,* rot u ub = 0, ie. 
ôu Pl/ôz — u P/âx = 0, (2.8) 
while*Y satisfies the two-dimensional wave equation 
(02/0x? + 02/02? — ve 202/o ES =À (2.9) 


together with the auxiliary relation div SY) = 0, te. 


Ou “dx + ôu_SV/ôz = 0 (2.10) 


The calculation is further restricted to the case where the force 
at the boundary is applied normally, then we have 


„= -fWòlx) 
where f(t) denotes the amplitude of the applied pressure 
(f(t) = O when —co <t <0) and ö(x) denotes the one-dimen- 
sional delta distribution. In view of the symmetry of the con- 


figuration with respect to the plane x = 0, u, is an odd Function 
of x, while u, is an even function of x, i.e. 


z,t)and u (x,z,f) = u_{—x, z, f) (2.12) 


limo T,, = Oand lim, ‚rt, (2.11) 


ulxzt)= —ul—x, 
As a consequence, u, and w_ need only be determined in the 
quarter-space ÛSx<oo, —oo0<y<oo, O<z<oo. Finally, it 
has to be taken into account that the waves generated by the 
source travel away from it. 


3. Method of solution 


We follow the steps as indicated in [2] and subject all time- 
dependent quantities to a one-sided laplace transform with 
respect to time. For example: 


Fs) = [5 expl—si)f(o)dr (3.1) 


in which s is a real, positive number large enough to ensure 
the convergence of integrals of the type (3.1). (It is assumed 
that the behaviour of f(f) and @(x, z, t) as f— oo is such that 
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an appropriate value of s can be found.) Let U =— Ux, z: s) 
denote the one-sided laplace transform of d = ü(x, z, f) and let. 


similarly, 7, etc. be the one-sided laplace transforms of t, 
ete. Then UP satisfies the differential equation 

(8°/dx* + d*/dz® — sfvp DUP 0 (3.2) 
together with the auxiliary condition 

U lgz — EU Slex =0 (3.3) 
while USV satisfies the differential equation 

(a /âx? + 0/ôz* — s°/v,®) US” — 0 (3.4) 


together with the auxiliary condition 
OU SYâx + AU Vaz = 0 (3.5) 


Further the constitutive relations (2.2) and (2.3) are replaced 
by corresponding relations involving the laplace-transformed 
quantities, whereas the boundary conditions (2.11) are re- 
placed by 
limo Ess 


= —F(s)ö(x) (3.6) 


Next we introduce the integral representations 


UP — (F(s)/2ri) Í (pi, + Yol )A(p) x 


_ 
EE 


= 0 and lim, , 7 


x exp —s(px + 7p 2)ldp (3.7) 
and 


Üen roa | il, — PI)BP) « 
x_expl—s(px + 7a2)ldp (3.8) 
in which 


Ves = (l/vp s° — p°)* with Re(y ps)20 (3.9) 


The right-hand side of (3.7) satisfies (3.2) and (3.3). the rieht- 
hand side of (3.8) satisfies (3.4) and (3.5); the choice of the 
square root in (3.9) will lead to waves travelling away from the 
origin. Substituting (3.7) and (3.8) in the boundary conditions 
(3.6) and using the representation 
Ò(x) = Ge | expl— spx)dp (3.10) 
iet 

we obtain two linear, algebraic equations for A(p) and B(p). 
The solution of these equations is obtained as 


Alp) = (1/28? — p°y/2uR(p); Bp) = py p/24R(p) GL) 
in which 
R(p) = (l/2ws* — p°)? + Pp? yes (3.12) 


At this point we observe that A(p) and B(p) have the following 
singularities in the complex p-plane: (a) branch points at 
p= £ l/r,, (b) branch pointsat p = + 1/vs, (c) simple poles 
atp= £ I/v, wherep = I/vg andp = —l/vg (vg = Rayleigh 
wave velocity) are the simple zeros of the right-hand side of 
(3.12). 


Now that UP and US” have been determined, the trans- 
formation back to the time domain follows by changing the 
path of integration in (3.7) and (3.8) to a curve along which 


(3.13) 


Since this process requires the integrands to be single-valued, 
branch cuts are introduced in accordance with (3.9), i.e. along 
Rel) = 0. Further, neither singularities nor branch cuts 
may be passed. The modified path of integration in (3.7) is 
then obtained asp = wolx, z, 1), together with its image p = w‚* 


Re(px + psa) =tT=>0 and Im(px + yps2) = 0 
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A 


with respect to the real axis, where w, is given by 


wp = (/r) sin(0) + i(t?/r2 — 1/vp°) * cos (0) (3.14) 


withT, <1 <oo and Tp = r/Vp 


The modified path of integration in (3.8) in general consists of 
the two parts p = wps(x, z, t) and p = wslx, z, v), together with 
their images p = @ps* and p = w;* with respect to the real axis, 
where ops Is given by 


Wps = (q/r) sin(0) — (1 lvs — tÀ/r?)* cos(B) + i0 (3.15) 
with Tos <t < Tand 
Tos = Tp sin(0) + Tl — vs-/rp°)* cos(B), 
T= rfvs 

and wm; by 

ms = (e/r) sin(0) + it?/r2 — 1/vg*)* cos(0) (3.16) 


with T, <t <00 


The part p = @pg, together with its image p = wps' with 
respect to the real axis, is only present for points of observation 
in the domain v/v, < |sin(0)| S l and is due to the presence of 
the branch points p = + I/v, (Fig. 2). Introducing in (3.7) and 
(3.8) 1 as variable of integration and taking into account the 
symmetry of the paths of integration with respect to the real 
axis, we obtain expressions for UP) and US of the form 


UP — FSG (x,z;s)and Us” — F(s) GEV (xx. zy s) (3.17) 


where 
GG = | exp(-— srt) g°'(x,2, 7) dr (3.18) 
Tp 
and 
Ïg 
Ge” 2 [ exp an sy) (x. z, r) dr d 
Tos 
| exp (— sr) 2“ (x,z, 1) dr (3.19) 
Ts 
with 
g© = (1/2) Im (pi, + 7p (Wp) À) A (Wo) (Oo p/ô7)] (3.20) 
g° = (l/r) Im [Gs @ps) ie Ops 1) BW) (Cops/ér)] (3.21) 


Im (p) 


Fig. A. Wave fronts of the P-wave, the PS-wave and the S-wave. 
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4 = (1/7) Im (ye (We) is) B (wo) (os/07)] 


(3.22) 


Expressions for u P) (x, z, ft) and SV (x, z, ft) are now obtained 
by applying the convolution theorem to (3.17) and using the 
uniqueness of the inverse laplace transform. The results are 
presented in Table | and Table 2. 


Table 1. P-wave as a function of time. 


Úiek Oz z‚t) 


VEA) g(x, z, 7D) dr 
P 


Table 2. SV-wave as a function of time (the term 4'** is only 
present in the domain ve/vp < |'sin(0)| < 1). 


sE f-DZ "Sz, de 
Ps 
ir es fC-DZPD x, 2, 0) de + 


+IEfA-Dg 2,0) dr 
5 


4. Discussion of the results 


The total elastodynamic wave motion @ generated by the line 
source is given byu =— úP + @*V, where u” follows from 
Table 1 and #'*Y follows from Table 2 (see Fig. 3). The P-wave 
is a cylindrical wave arriving at : = T‚, where T, is the travel 
time for P-waves travelling along a straight line from the line 
source to the point of observation. The SV-wave consists of a 
PS-conversion wave arriving at f = Tpg and a cylindrical SV- 
wave arriving at f = T. The PS-conversion wave is only present 
in the domain va/vp < | sin(0)| £ l and precedes the cylindrical 
SV-wave. Tpg is the travel time for SV-waves travelling from 
the line source a distance x—z (wa/vp) (l—vs°/vp°) “along the 
boundary with P-wave velocity and then along a straight line 
to the point of observation with S-wave velocity; Ts is the travel 
time for S-waves travelling along a straight line from the line 
source to the point of observation. At the boundary of the 
solid there is a singularity in the displacement vector travelling 
with velocity vj; this singularity is due to the presence of the 
poleatp = + I/vp. 

The features that are characteristic for the configuration 
occur in gg PS and g°. The numerical computation of 
these quantities involves only algebraic operations of the type 
complex. 
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Transiënte verschijnselen bij synchrone machines 


door prof. dr. ir. F. A. W. van den Burg, Afdeling der Elektrotechniek, T.H. Delft 


Summary: Ön transients of the svnchronous machine. 


Frequently, transients of the synchronous machine have been studied by analytical methods and pictured almost 


exclusively by the phase portraits of the load angle. 


A different portraiture is introduced with attractive features as to the insight in the dynamical behaviour due to 


impact loads. 


Besides, the validity of the second order differential equation, which is called “the differential equation of the 


synchronous machine’, is discussed. 


|. Inleiding 


De afbeeldingen die het stationaire gedrag van elektrische 
machines beschrijven behoren tot de standaardafbeeldingen van 
de vakliteratuur. Tot zulk soort afbeeldingen kunnen voor de 
synchrone machines de koppel-lasthoekkarakteristieken en de 
stroomdiagrammen worden gerekend. 

Stelt men voor de synchrone machine, onder niet te essentieel 
beperkende voorwaarden, het volledige stelsel op van differen- 
tiaalvergelijkingen dat het elektrodynamische gedrag van de 
machine geheel beschrijft, dan vindt men een stelsel van de 
vierde orde voor het geval dat er geen demperkooi aanwezig 
is en een van de zesde orde indien er wel een demperkooi 
aanwezig is. Na een passende assentransformatie kan men deze 
stelsels als autonome stelsels schrijven. Dan kunnen ook de 
stationaire oplossingen van de oorspronkelijke stelsels van de 
vierde of de zesde orde worden gevonden en de reeds genoemde 
karakteristieken en diagrammen worden geconstrueerd. 

Voor het verkrijgen van afbeeldingen van overgangsver- 
schijnselen die optreden na plotselinge veranderingen die op de 
machine worden uitgeoefend, kan men in principe twee geheel 
verschillende wegen inslaan. Men kan de genoemde stelsels 
direct bewerken met behulp van rekentuigen of men kan deze 
stelsels zodanig analytisch bewerken dat enig inzicht op grond 
van de verkregen mathematische vereenvoudigingen mogelijk 
wordt. 

Tussen de beide uiterste wegen ligt de middenweg van de ge- 
doseerde analytische voorbewerking van het probleem. Voor een 
daarop volgende bewerking met een digitaal rekentuig betekent 
zo'n voorbewerking veelal het traceren van de sporen waarlangs 
de programmering zich voor een economisch gebruik van dat 
hulpmiddel moet bewegen. Zo'n middenweg leidt bij het gebruik 
van analogonrekentuigen bijna steeds tot aanzienlijke compli- 
caties voor de daarbij vereiste programmering. 

De eerst in de afgelopen decennia ingezette snelle ontwikkeling 
van beide rekentechnische hulpmiddelen maakt het begrijpelijk 
dat vóór 1940 vooral de analytische bewerking tot op het niveau 
van een mathematisch interpreteerbare formulering de volle 
aandacht had. Reeds in 1911 leidde Dreyfus [1] een benaderende 
differentiaalvergelijking van de tweede orde af voor de (twee- 
fasige) synchrone machine. Aan zijn afleiding werden vele 
details toegevoegd, bijvoorbeeld door Bödefeld [2] Edgerton en 
Zak [3], Edgerton en Fourmarier [4]. Lyon en Edgerton, [5]. 
Timascheff [6]. Hun afleidingen berusten op gevoel voor ver- 
houdingen, fysisch inzicht, intuitie, of hoe men dit ook noemen 
wil, maar zeker op een grondige praktische ervaring met deze 
machines. Toch moet worden gesteld dat de voorwaarden waar- 
onder de afgeleide differentiaalvergelijkingen met een te schatten 
nauwkeurigheid geldig zijn, volledig in het duister bleven. Mede 
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daarom werd reeds in 1935 het gehele stelsel differentiaalver- 
gelijkingen met behulp van een analogonrekentuig bewerkt [7]. 
‚ Een goede revisie van de mathematische formuleringen van de 
genoemde stelsels brachten, enigszins onafhankelijk van elkaar, 
Laible [8] en Concordia [9]. Op mathematisch gebied trok daar- 
naast de toepassing van de theorieën van de niet-lineaire 
differentiaalvergelijkingen op de (door Dreyfus afgeleide) ver- 
gelijking van de synchrone machine via Vlasov [10], Minorsky 
[1] en Stoker [12] zo sterk de aandacht dat het soms tot nieuwe 
mathematische beschouwingen aanleiding gaf [13], maar zeker 
bewerkte dat de afbeeldingen van de overgangsverschijnselen 
die bij de synchrone machine optreden, vooral door middel van 
een faseportret van de lasthoek werden weergegeven. 

Dat andere afbeeldingen vaak een beter inzicht in het gedrag 
van overgangsverschijnselen geven zal worden geïllustreerd door, 
voor een viertal stootveranderingen op de synchrone machine, 
naast de gebruikelijke faseportretten van de lasthoek, een on- 
gebruikelijk type van afbeelding van hetzelfde viertal verande- 
ringen weer te geven. 

Deze afbeeldingen zijn met behulp van een analogonrekentuig 
verkregen. Voor de met het gebruik van dat hulpmiddel weinig 
vertrouwde lezer zij gewezen op de ‘Vakantieleergang Elektro- 
techniek 1959 [14] en op de recente dissertatie van Köppel [15]. 

Tegen de bovenstaande achtergrond zal eveneens, uitgaande 
van een formulering overeenkomstig aan die van Laible en 
Concordia, nog een analytische voorbewerking van het be- 
schrijvende stelsel van differentiaalvergelijkingen worden be- 
handeld waarmede de door Dreyfus reeds afgeleide benaderende 
beschrijving direct wordt gevonden. 


2. Vergelijkingen 


Voor de synchrone machine met gladde rotor zonder demperkooi 
(zie fig. 1) waarvan de stator is aangesloten op een star symme- 
trisch driefasennet vindt men als beschrijvende vergelijkingen 
in een genormeerd dq-assenstelsel dat op details na overeenkomt 
met het door Laible [8] gebruikte assenstelsel, 


E = PJ+(1-ò)GP+b:D=Klit,=j,d+JGP (I) 


K= [Xpo A =EI Dix = ej, @) (2) 

E = [—sin (Ö), + cos (ò), + e‚] (3) 
1 1 1. k 

P= ë Geit - EES Hel 
at, de en | bk «1 

a= rik, a, = FIX, (5) 
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Fig. 1. Synchrone machine (motor); tekenafspraken. 


r __— statorweerstand | (6) 
X__= synchrone reactantie (7) 
r‚ = weerstand veldwikkeling (8) 
X, = reactantie veldwikkeling bij netfrequentie w,/2r (9) 
Ö — lasthoek,ò =dö/dr,t =wt,t=tijdinsecunden (10) 
k _—= koppelingscoëfficiënt, « = 1 —k* = lekcoëfficiënt (11) 


E. resp. J, resp. ® = getransformeerde en genormeerde excitatie- 
spanningsmatrix, resp. stroom-, resp. fluxmatrix, (12) 


! 


Arm 
In (1) en (2) betekent X’ de matrix die uit X door transpositie 


ontstaat. Voorts zullen nog worden gebruikt de met « en a, 
samenhangende hoeken Z en Z, volgens, 


cos (5) = X/Z, sin (&) =r/Z,Z = (r° + X°)*,a = tan (5) (15) 


genormeerd asmoment (13) 


genormeerd polair traagheidsmoment van de rotor (14) 


cos(4,) = X,/Zy, sin (&‚) =r,/Z,,Z,= (‚+ X,°) ê 

dr = {an (&,) (16) 
Enig idee van de orde van grootte van de gebruikte constanten 
geeft tabel 1. | 


Tabel 1. Gebied van de machineconstanten 5, a, #,, j,. 


Machine 


_ Turbo- 0,020 0,013 1500 


generator | 0,030 0,035 


0,010 0,035 150 


0,030 0,150 500 


De vergelijkingen voor een machine met uitgebouwde polen en 
met een volledige demperwikkeling verschillen niet principieel 
van (1). Zo moet men dan bijvoorbeeld in de plaats van (3) en (4) 
stellen: 
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Ki 


d, 
p=. Dn (17) 
JN 
at, | 
| Kin ki 
l ko 
K= | Kap | kor 
ko l 
ks Kos | 


Voor het verdere betoog zijn deze verschillen niet van wezenlijk 
belang. 


3, Stroomdiagram 


Het stelsel (1) laat een particuliere oplossing toe met constante 
x= x°, (x = ej, p, ò), die samenhangen volgens 


E= PJ 4+Gb; D=KJI; (== JGD =kijfjf (18) 


Met (17) kan men voor deze stationaire toestand allerhande 
van technisch belang zijnde grootheden uitdrukken in ò°, de 
constante lasthoek en in j,°, de constante bekrachtigingsstroom, 
waarbij in deze genormeerde grootheden de combinatie kj, 
voorstelt dat deel van de magnetomotorische kracht van de 
veldwikkeling dat maximaal met een d- of q-statorfase kan zijn 
gekoppeld. 

De onderlinge samenhang van de bedoelde grootheden wordt 
weergegeven in het bekende cirkeldiagram van de synchrone 
machine (zie fig. 2). Dit diagram zal tevens worden gebruikt 
voor het afbeelden van overgangsverschijnselen. Daarom 
worden nog gedefinieerd : 


Es (e‚ +e,°)°, (n.b:e,=e, = Dij,= Os +) (19) 
PS =eada Hed, = —Îasin (0) + j,cos (0) = j, cos (0) (20) 
(21) 


waarmee, juist voor niet-stationaire toestanden, het aan de 
stator toegevoerde wattvermogen p‚ en het blindvermogen 
p®* als functie van de tijd r zijn gedefinieerd en ook de fase- 
hoek p = @ (rt) is vastgelegd. Deze definities geven de beschik- 
king over paren van variabele grootheden die in de stationaire 
toestand als constante grootheden bij elkaar behoren en van 
techmisch belang zijn. 

Het diagram van fig. 2 verschilt slechts van het gebruikelijke 
stroomdiagram door het gebruik van genormeerde grootheden 
en door een draaiing van het gebruikelijke diagram over 
+ n/2en wel + of — afhankelijk van het gekozen standaard- 
werk waaruit men zo’n stroomdiagram ontleend kan denken. 
Ook het aflezen van het (constant.) asvermogen p‚; en van de in 
de stator gedissipeerde warmtestroom p‚“"“ komt in fig. 2 op een 
niet gebruikelijke manier tot stand. In de gekozen stand kan 
men tenslotte, door de aard van de normering, het diagram op- 
vatten als een stroom-, een vermogens- of als een admittantie- 
diagram. 


pe = Erle Edda —j4 cos (ò) —Ja sin (ó) = j, sin (p) 


_=i ZO ak 5 
| Ege generator netwerk motor =/ hd 
0 of B, ee e 
/ RR 7 Ei 5: 
rs B" f B nf Cn, / à susceptantie 
id (A ke 5 j br 
rr Of be Se c' 
EF w Of / Ta, / 0 R 
| ĳ f KM „ _ conductantie 
/ ZN EN Mn 
erk oeh ramt mese: de capacitief 
/ 4 | l € Ah wattloos- verm. 
/ | a hr A Dr 
ii Lef TOA ER 
| Es} ) / 0 | SON Oog? 
- id an = B Dn Ps 
/ MeT C \ watt-verm. 
‚_\N | 
Nd 
| ÁÀ 
| È | [o Jeostai A sini 
nn. \ os) FA N=cos{erl | sin); -coslä) 
N \ 4 157 
Nen 
À ie 0 ZEN 
S 


Er Vn 


An 


Fig. 2. De synchrone machine; het genormeerde cirkeldiagram voor 
de stationaire toestanden, 


4. Afbeelding van stootveranderingen 


Uit (1) tot en met (4) volgen voor 


—Ô4 = Hg — (1-0) 0, + sin (Ö) (22) 
—Ô,= 4, + (1-0) p, — cos (ò) (23) 
Op dy Et =O (ijd) (24) 

Pa = da + hip Pe, = In Pf = ha + Jo (25) 
-j,ò= kij, t8 (26) 


de vergelijkingen voor de machine met een constante bekrach- 
tiging en een constant lastmoment. Het stelsel (19) tot en met 
(26) werd, voor een 3-kVA motor met 


u = 0,03; a, = 0,035; a = 0,30; j, = 250 (27) 


met behulp van een analogonrekentuig gesimuleerd. 

Uitgaande van de begintoestand (p‚; = t° = 1,44; kj,“cos(a) 
= 2) werden de dynamische karakteristieken geregistreerd die 
het verband aangeven tussen enerzijds p‚’'(t) en p‚(t) (of, met 
andere woorden, tussen j‚(t) en (rt), zie fig. 3) en anderzijds 
tussen kj, (t) cos (2) en ò (t) (zie fig. 4), en wel na het aanbrengen 
van de laststoten (geval a en b), alsook na het aanbrengen van de 
bekrachtigingsstoten (geval c en d). 

Omdat het stelsel (19) tot en met (26) een differentiaalver- 
gelijking van de vierde orde omvat, kan niet worden verwacht 
dat de overeenkomstige dynamische karakteristieken van de 
fig. 3en 4 elkaar dekken. Toch blijkt dit in hoge mate het geval 
te zijn. Dit betekent dat het genoemde stelsel redelijk moet 
kunnen worden benaderd door één (niet-lineaire) differentiaal- 
vergelijking van de tweede orde. 

In fig. 5 zijn de bij de gevallen a, b, c en d behorende fase- 
portretten van de lasthoek ò weergegeven. De onderscheiding 
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van het subtransiënte gebied — waarin ò bij benadering constant 
blijft — en het transiënte gebied is voor geval c en d uit de fig. 3 
en 4 nog mogelijk, maar niet in de fig. 5. De afbeeldingen die zijn 
gesuperponeerd over het cirkeldiagram geven dan ook meer 
informatie en een beter inzicht in het gedrag van de overgangs- 
verschijnselen. 


5. Benadering 
Uit (1) volgt voor E: 

d | 
E= P+GK)I OG —e)K]J;e =diag[1, 1, 1] (28) 


De schrijfwijze van (28) is zodanig ingericht dat voor de statio- 
naire toestand met ò — ò“ en (o.a) J = J' de tweede term uit 
het rechter lid van (28) verdwijnt. 

Men noemt J'°' uit 


E= (P+ GK)J” (29) 


de quasistationaire beschrijving van de spanningsvergelijking 
(28) en het daarmede bepaalde moment 


Le, =j,d + kj} COS (@) [sin (Ö + 2) —kj, sin (2)] (30) 
de quasistationaire momentvergelijking. Als e f = e‚ en 
ts = fo, kan men (30) vervangen door 

j,ò + kj,“ cos (5) [sin (+2) — sin (ö° + | =O GĲ) 


een vergelijking die de ongedempte beweging van de machine 
beschrijft en aanleiding geeft tot faseportretten van ò + # die 
door Stoker [12] uitvoerig zijn behandeld. 

Men kan daarom verwachten dat als men voor J in de tweede 
term van (28) de J'” uit (29) substitueert, de dan te bepalen 
stroommatrix uit (28) een dempingsterm in de momentverge- 
lijking introduceert. 

Stelt men nu 


A = (P+GHO!; D= GG) K (32) 
T 


mm 


| 
| 
| 
4 


ú 


| 
| 
| 
| 


en :! 
Fig. 3. Dynamische karakteristieken van (p‚'; p°), dus ook van 
(j,: @), uitgaande van [p‚S = 1,44; kij cos(%) = 2], en wel: 

a. met sprong van p‚‚ van 1,44 naar 0,96; 

b. met sprong van p‚; van 1,44 naar 1,94; 

c. met sprong van kj; cos(Z) van 2 naar 15; 

d. met sprong van kis cos(Z) van 2 naar dn 
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A 


+3 


Fig. 4. Dynamische karakteristieken van (kj, cos(?); ò), uitgaande 
van [p‚S = 1,44; kjf cos(Z) = 2], en wel de gevallen a,b, c en d (zie 
fig. 3). 


Fig. 5. Fasediagrammen voor de onder fig. 3 en 4 vermelde gevallen 
a, b, cen d. 


dan is A een constante matrix en D een matrixoperator en kan 
men (28) schrijven als 


J= AE + ADJ (33) 


Door herhaalde substitutie van het rechterlid van (33) voor J 
uit dat rechterlid vindt men zo de ontwikkeling 


Je > (AD) AE (34) 
n= 


waarvan de partiële som 
‚_N 


JOY (AD)” AE (35) 


n=Ü 
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voor N = 0 de quasistationaire J'® uit (29) levert. Voor N=l 
vindt men tenslotte als benadering van de momentvergelijking: 


1 | 
Ji ò + m k’ cos° (@)[sin° (Ô +&) — 2, «sin (ô +#)] oe 
XK 
36) 
(kj) re ô+ kj} cos(@)[sin(ô +#) —sin (6 +a)}=0 
een vergelijking die ook door Bödefeld [2] is afgeleid en waarvan 


de faseportretten van ò + & in Minorsky en Stoker uitvoerig 
worden behandeld. 


6. Een bovengrens voor | ö | bij het toepassen van (36) 


Met de notaties 


==; X= w(dxjdö) = wr; (x= ò,Ô,ô,…) (37) 


kan men (34) met behulp van (32) en (37) schrijven als, 
En d n} 

J=Y [W(AGK- AK a) AE (38) 
n= 


Men kan nu, de methode uit Bellman [16] volgende, elke com- 
ponentenvergelijking van (38) majoreren door één scalaire 
reeks; daarbij wordt gebruik gemaakt van het begrip ‘norm’ 
van een matrix B, genoteerd als || Bj]. Gebruikt zal worden de 
euclidische vorm van zo’n matrix met als definitie 


|B} = [som van de kwadraten van alle elementen van B] * (39) 
Bepaald werden nu 


| 
|All= En [I + (k° +245) cos° (5)}° =a, = 22 xr, (40) 


fi 
ES 2 2 
AK +42 0 4 2E 02) eos? (@) | = 
% | | 


= A, 2e (41) 


55 H 
LAGK || = cos (£) 2 +k? cos° | =a, 3? (42) 
waarbij voor de benaderingen van a, a, en a, de gegevens uit 
tabel 1 werden gebruikt. Dan kan men volgens de genoemde 
methode naar aanleiding van (38) stellen 


af Cn) 
i,lsa y lol(as Haden IEI, e=daf) 43 
n= 


waarbij de absoluuttekens om d/dò aangeven dat bij de diffe- 
rentie naar ò van Een van w alleen de absolute waarden van de 
resultaten tellen. 


In het bijzonder moet men in (43) lezen 
Î 


} F 


d cm de, - de \ - de A3 2 
| IEl=|{ sl) + Ù + :) (m4 0) 
dô | dó” dA” dô” + 
(44) 
d 5 a 
eri lEl=| e+ ele | (45) 


ti |= 


hetgeen voor E = [—sin (Ö), + cos (Ö), + e‚°J geeft: 
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| ed 


d Kar 


L—| (voor m = 0) 
dò 


WE|=1 voor m +Oen =|. He 


(46) 


Opdat het rechterlid van (44) beperkt blijve en daarmede de 
schatting volgens (44) aan een eindige bovengrens gebonden 
blijft, moet elke bijdrage aan (44) beperkt blijven en dus ook 
de bijdrage die ontstaat door alle termen die @', »'', 1''’ ES 
bevatten, te verwaarlozen. Door de bedoelde bijdrage af te 
zonderen kan men voor (44) schrijven, 


AES:D) ola, val Sj” Diens 5 
n= Ö 


verdwijnende functie. 


Opgemerkt wordt dat het gedrag van @', @'', @''’, …. sterk 
af zal hangen van de grootte van het polaire traagheidsmoment 
J, en het is niet onredelijk te onderstellen dat de eerste term 
van (47) voor grote j, de belangrijkste bijdrage zal vormen 
voor de maximale excursies van ». Nu geldt 


an 


ie dt en „l,d © 
bà Ko ei eerd IEll=) |a; wi 5 + 


n= n= 


a dw | a,\ dw 
+ JIE Ee )(e el (IEl= 
a/ dô a, / dô 


ee En ES 
=d la as +(+ €”)? 1) (48) 


n= 3 
een bijdrage die slechts eindig blijft voor 


I 
3422 xy 


Voor ax, = 0,15 (zie tabel 1) zal dus (34) alleen voor Löl < 0,05 
kunnen convergeren. 


|Öl=|w|< 


et IA DAsE 


550.12 


A note on unified field theory 


bij prof. dr. R. Timman, pet University of Technology, Department of General Sciences, 


Subdepartment of Mathematics 


Summary: Á unified field theory is considered, where the electromagnetic field is represented by the anti-symmetric 


part of a non-symmetric connexion. 


Most attempts at a unified field theory (e.g. V. Hlavaty [1]) 
introduce a non-symmetric metric tensor and relate the gravita- 
tion to the symmetrical part and the electromagnetic properties 
to the antisymmetrical part. The first gives an ordinary Riemann 
space and its tensor is usually used to form covariant vectors. 
Apparently both tensors are, right from the beginning, com- 
pletely separated. The existence of electrodynamics in special 
relativity, where the space-time continuum has a constant 
metric tensor, while the electromagnetic field is governed by 
the antisymmetric field tensor leads to a concept, where the 
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De moeite die men aan analytische arbeid moet besteden om 
benaderde vergelijkingen voor het stelsel (1) te verkrijgen is, 
ten opzichte van het gemak waarmede men met behulp van een 
analogonrekentuig een goed inzicht in het gedrag van het niet 
benaderde stelsel verkrijgt, onevenredig groot en geeft boven- 
dien resultaten waarvan de nauwkeurigheid nauwelijks valt te 
schatten. 
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metric tensor is symmetrical, but the connexion ’‚ has an 
antisymmetric part 


Sj E î Fi zi Ti) 


By inspection of the transformation laws for the connexion 
te.g. Schouten [2) this object appears to be a tensor. which is 
called the torsion of the space. 

Here we shall derive a Hamilton principle in this space, 
which allows us to postulate a relation between the torsion 
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tensor and the electromagnetic tensor. A similar consideration 
in completely different notation is given by Sciama [3]. 


The parallel displacement of a vector vi, which is transferred 


from a point x‘ to a point x* + dx‘ in the space with holonomous 
coordinates x' is defined as a displacement for which 


òv' =de’ + Fi, vt dx? = 0 


The object Ti is called the connexion in the space and the 
vector ôv' is called the covariant differential of v°. It is a linear 
homogeneous function of the displacement dx’ and it satisfies 
the ordinary rules with respect to addition and multiplication 
of a differential. As was mentioned above 15, need not be 
symmetrical and 


Sin = 3 Pie Vi) 

is a tensor, the significance of which we shall investigate. If in 
a point x! we have two line-elements dx’ and òx' and we dis- 
place òx! along the line-element dx’, we have for a parallel 
displacement the endpoint 

xt + dei + Ox! + Tij xt dai = x° + dx! + òx! + dôx' 

on the other hand for a parallel displacement of the vector dx! 
along the vector òx/ we have 


xi + ôx! + dx! + Di, ôx! dri = xi + öxl + dx! + ödx' 


In our case the endpoints do not coincide, but they differ by a 


vector 


döx' — âdx' = (Ti —_ Fi j)öx' dx! = 2 Si öxk da). 


This property will be used for the derivation of a Hamilton 
principle. 
H 


| ds = stationary 


A 


for different curves, running from a fixed point A in space-time 
to a fixed point B. Here ds is the line-element, defined by the 
symmetrie metric tensor g,, 


ds =de de” 


It can be shown (Schouten [p. 132, 2}), that, if we put the co- 
variant derivatives 


VI = 0 
that the tensors g,, and S, define uniquely the connexion. 


Fien} Si — Sijn +Sijn 


where the tensors S;, and S, are derived from S;,’ by raising 
and lowering of indices 


ik I Iam Sit 5 Sie SG Gum Sij 


and |; |} stands for the Christoffel symbol 


Lits 


En 
ijkf — 24 


rd e Ì 
Oda + OI — CA jks 


Consider now in this space the variation of an integral 


B 


Í L(x', u) ds 


„A 


‚ dx 

where U = — 
ds 

a variation òx'(s) with a corresponding àu' — 8 


‚ In each point x'(s) of the extremal we define 
dx 
ds 
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From the relation derived above, we can infer that, remarking, 


thats means a differentiation ge the line-element dx’ 


| de d x 
duld =S 8x 
ds ds 5 5 
The condition for stationarity of the integral now reads 


BH 


B 
| öL ds | IK (xl 4 Bel, ul + Bu) L(x', 1} ds = 
4 A 


de 8x + ee du’) ds =0 
Öx' Öu' 
Ee | 


for arbitrary variations òx' which vanish at the endpoints. 
Substitution of the expression for òu' gives 
B 


ÔE. , OE Hd. > êL Te 
EK | òx u’ \ds=0 
âx! ôu' ds du E 


A 


or after partial integration of the second term, 

r/êL d êl 

(ESE Bejven 
dx ds êu' êu' 

A 


d êL Ll ôL 
d- 


ds êu' ôx' êu’ 


We now remark, that L = /9;; w u“, which, along the extremal 
has the value 1 and carry out the differentiation. The first two 
terms give, similar results as those for a symmetrical connection, 
but the last term gives an addition. 


ds 


Î X 
Bi + De | uu — 25, A u ul =0 


From this result we find the equation for the geodetic line 


For the equation of motion of an electrical particle in a gravita- 
tional field we have (Landau and Lifschitz [4]) 


du’ HE Rak l ok 
Stel u u me Fe ü 


where F, is the electromagnetic field tensor. 
This leads to a tentative identification 


ij 
BE F‚. ze Sj ä, 
mc 
Regarding the formula for the vector potential A, we have 
FV }A,- VrAj= 0, A, — OA, TA, + Dij= 
—2Si,A, 


This shows, that in the case. of a flow field, which satisfies the 
relation 


Al nt —= 
Öju, — Ô,u,= 0 
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we may indentify u, with the vector potential 
pr} 
mm 
A, En Re U, 
! 
It might be of some interest to work out further consequences 
of this hypothesis. 
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